S645.0998 



The conformal window in QCD and supersymmetric QCD 



Einan Gardi and Georges Grunberg 



Centre de Physique Theorique de l'Ecole Polytechnique* 
91128 Palaiseau Cedex, France 
email: gardi@cpht.polytechnique.fr, grunberg@cpht.polytechnique.fr 



Abstract 

In both QCD and supersymmetric QCD (SQCD) with Nf flavors there are conformal 
windows where the theory is asymptotically free in the ultraviolet while the infrared physics 
is governed by a non-trivial fixed-point. In SQCD, the lower Nf boundary of the conformal 
window, below which the theory is confining is well understood thanks to duality. In QCD 
there is just a sufficient condition for confinement based on super convergence. Studying 
the Banks-Zaks expansion and analyzing the conditions for the perturbative coupling to 
have a causal analyticity structure, it is shown that the infrared fixed-point in QCD is 
perturbative in the entire conformal window. This finding suggests that there can be no 
analog of duality in QCD. On the other hand in SQCD the infrared region is found to be 
strongly coupled in the lower part of the conformal window, in agreement with duality. 
Nevertheless, we show that it is possible to interpolate between the Banks-Zaks expansions 
in the electric and magnetic theories, for quantities that can be calculated perturbatively in 
both. This interpolation is explicitly demonstrated for the critical exponent that controls 
the rate at which a generic physical quantity approaches the fixed-point. 
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1 Introduction 



In multi-flavor QCD, there is a conformal window [1], namely a region of Nf values 
for which the theory is asymptotically free at short distances while the long distance 
physics is governed by a non-trivial fixed-point. This is a non-Abelian Coulomb 
phase in which quarks and gluons are not confined. 

The upper boundary of the conformal window is determined according to the 
sign of the (3 function 

PW = d^ = -(/W + /W + -"), (1) 

at small coupling x = a s /7i. When the first coefficient of the perturbative f3 function 
[2], 

A = i(^-^), (2) 

changes its sign the theory changes its nature from the asymptotically free phase 
r = Nf/N c < 11/2, to the infrared free phase r > 11/2. The transition point in 
N c — 3 QCD is at Nf = 16.5. For r < 11/2 the f3 function is negative for a 
vanishingly small coupling (/3 > 0), but due to the second term which has an 
opposite sign < 0), /3(x) reaches a non-trivial zero at some x FP ~ —/3o/j3± > 
[3]. x FP approaches zero as r approaches 11/2, and then quarks and gluons are 
weakly coupled at all scales. Finally, the smallness of x PP justifies the use of the 
2-loop (3 function. 

On the other hand, the lower boundary of the conformal window, below which 
confinement sets in, is much harder to tackle. One approach to confinement, the 
so-called metric confinement [4], defines confinement as a phase in which transverse 
gauge field excitations are excluded from the space of physical states which is defined 
through the BRST algebra. It was shown in [4] that as long as a certain condition is 
obeyed by the gauge field propagator, metric confinement is implied. This condition 
is most conveniently expressed in the Landau gauge, namely that the absorptive 
part of the gluon propagator D(Q 2 , fi 2 , g) (in this gauge) is superconvergent, 

/ dk 2 p(k 2 ,p 2 ,g) = 0, (3) 
No- 
where fj 2 is the renormalization scale and p(k 2 ,p 2 ,g) = (l/7r) Im {D(— k 2 , p 2 , g)}. 

Assuming analyticity of the gluon propagator D(Q 2 , p 2 , g), the superconvergence 
relation (3) was shown to be a direct consequence of renormalization group invari- 
ance, provided D(Q 2 , p 2 , g) vanishes faster than 1/Q 2 at large Q 2 . Due to asymp- 
totic freedom, the last condition depends only on the sign of the 1-loop anomalous 
dimension of the propagator (in the Landau gauge), given by 

1 /13 2 \ 
7oo = — A [ Y N C - -N f ) . (4) 

If 7oo is negative the superconvergence relation (3) holds. 
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Thus, in this approach, a sufficient condition for confinement is that 700 < 
and therefore the lower boundary of the conformal window cannot be lower than 
r = 13/4. If super convergence is also a necessary condition (this has not been 
shown) then the phase transition should be at r = 13/4. In iV c = 3 QCD this 
corresponds to Nf = 9.75, i.e. between 9 and 10 flavors. 

A natural question to ask is whether the super convergence condition necessarily 
implies that also quarks are confined. To the knowledge of the authors no complete 
answer has yet been given to this question, although it has been shown that the 
super convergence criterion is consistent with a potential that is approximately linear 
in some intermediate scales [5]. 

There are several other approaches to study the phase structure of multi-flavor 
QCD, such as the instanton liquid model [7], the gap equation [6, 8] and computer 
simulations on the lattice [9, 10]. The new lattice results [10] are inconsistent with 
old ones [9] and with the superconvergence criterion for confinement: they indi- 
cate that the phase where no confinement nor chiral symmetry breakdown occurs, 
stretches down to Nf = 7, and thus only for 6 flavors and below QCD appears as 
a confining theory with chiral symmetry breaking, as we know it in the real world. 
In spite of these contradicting evidence, we assume in this paper that the bottom of 
the conformal window is as implied by superconvergence. 

Recently the presence of a fixed-point in QCD was studied as a function Nf by 
considering the perturbative (3 function [11, 12]. Three approaches where consid- 
ered: a direct investigation^ of the equation (3{x) = in physical renormalization 
schemes [13], the Banks-Zaks expansion [1, 14, 15] and the analyticity structure of 
the coupling constant. 

The direct investigation of zeros in the QCD (3 function in physical schemes 
[11] shows that at 3-loop, a fixed-point can appear for most effective charges above 
Nf ~ 5. However, presence of a fixed-point at the lower end is very sensitive to 
higher-loop corrections, and thus cannot be trusted. 

The Banks-Zaks expansion is an expansion in the number of flavors down from 
the point where (3q changes its sign. The expansion parameter is proportional to 
Po, or in N c = 3 QCD to (16.5 — Nf). It was found in [11] that Banks-Zaks series 
for different QCD observables behave differently: in some cases the coefficients are 
small and the expansion is reliable and in other cases it seems to breakdown at order 
(16.5 — Nf) 3 already around 10 or 12 flavors (see fig. 6 in [11]). Here we further 
interpret these results. 

In real-world QCD there are Landau singularities in the perturbative coupling, 
which signal the inapplicability of perturbation theory to describe the infrared 
physics. These non-physical singularities are usually assumed to be compensated by 
non-perturbative power-like terms in any physical quantity. Thus causality is real- 
ized only at the non-perturbative level. The existence of a perturbative fixed-point 
opens up the possibility that the perturbative coupling will have no Landau singu- 
larities [11, 12]. We will assume that within the conformal window non-perturbative 

'''The relevant refs. appear in [11]. 



2 



effects are not important as long as they are not implied by perturbation theory that 
is as long as the perturbative coupling is causal and small. Note that the causality 
requirement is stronger than the requirement to have no space-like Landau singu- 
larity. 

The simplest example where causality of the coupling can be achieved at the 
perturbative level, without additional power corrections, is the 2-loop coupling. In 
[12] an exact explicit formula for the 2-loop coupling as a function of the scale was 
introduced, which enabled a complete understanding of the singularity structure of 
the coupling in this approximation. It turns out that the condition for causality of 
the 2-loop coupling is c = Pi/Po < —fto- This condition translates in N c = 3 QCD 
to N f > 10. 

As stated above the lower boundary of the conformal window implied by super- 
convergence is also between 9 and 10 flavors [4]. Thus basing on the superconver- 
gence criterion, we find that the 2-loop perturbative coupling is causal in the entire 
conformal window. This suggests that the perturbative analysis in the infrared is 
reliable down to the bottom of the window. On the other hand, perturbation theory 
cannot describe the infrared physics in the confining phase. We therefore intend 
to study more carefully down to what Nf can we trust perturbation theory in the 
infrared, and in particular, when does perturbation theory signal its own inapplica- 
bility. The very same questions can be asked also in supersymmetric QCD (SQCD), 
where more is known about the phase structure. We therefore study here both QCD 
and SQCD and compare the two. 

A few years ago Seiberg lead a revolution in the understanding of supersymmetric 
gauge theories 1 ''. Of particular interest to us is the phase structure of N — 1 SQCD 
in which non-Abelian electric- magnetic duality plays a major role [16]. The general 
arguments available in the supersymmetric case do not apply in the absence of 
supersymmetry and in fact the phase structure of a supersymmetric theory can 
be quite different from that of its non-supersymmetric parallel [17, 18]. Still the 
comparison can be very enlightening. 

SQCD, just like QCD, has a conformal window where the theory is asymptoti- 
cally free, and is governed by a non-trivial fixed-point in the infrared. The picture 
described in [16] is the following: in the upper part of the conformal window the 
fixed-point value of the coupling is small, and thus the theory is weakly coupled at 
all scales. The massless fields that appear in the Lagrangian conveniently describe 
the physics at any scale. When Nf becomes smaller (for a fixed N c ) the theory 
becomes strongly coupled in the infrared. Then, it does not make sense anymore 
to describe the physics in terms of the original massless fields. Nevertheless, in 
the infrared limit the theory has an effective description in terms of a dual theory: 
starting with an original supersymmetric theory with an SU(N C ) gauge symmetry 
and Nf massless chiral quark superfields (Qi, % = l,2,...iVj) and their anti-fields 
(Qi, i = 1, 2, ...Nf) the dual theory is an SU(Nf — N c ) gauge theory with Nf chiral 
quarks superfields i = l,2,...Nf) and their anti-fields (q^ i = l,2,...Aj), and 

*For recent reviews see [17, 18]. 
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an additional superpotential describing a Yukawa interaction between color-singlet 
mesons and the quarks superfields: 

W = y/\M t j q i q>. (5) 

The relation between the theories is referred to as duality since the dual of the 
dual theory is again an SU(N C ) gauge theory. In the conformal window the dual 
theory, just like the original one, is asymptotically free and has a non-trivial infrared 
fixed-point. Contrary to the original theory, the dual theory becomes weakly coupled 
as Nf decreases, until the point where it becomes infrared free. Since the dual theory 
is weakly coupled when the original one is strongly coupled, and vise-versa, Seiberg 
refers to this duality as a non-Abelian generalization of the electric-magnetic duality. 
The duality picture is valid also outside the conformal window, where one of the 
theories is infrared free and the other is confining, but here we concentrate on the 
conformal window. 

In [19] (see also [20]) it was shown that the lower boundary of the conformal 
window implied by the superconvergence criterion for confinement in SQCD, coin- 
cides with the lower boundary implied by duality which is the point where the dual 
theory becomes infrared free. This gives additional support to the whole picture, 
with the advantage that the superconvergence criterion can be applied also in the 
non- sup er symmet r ic it was originally done in [4, 5]. 

According to Seiberg's description, in SQCD the electric theory is strongly cou- 
pled at the bottom of the window. We shall demonstrate that this strong coupling 
behavior manifests itself already at the perturbative level, through appearance of 
Landau singularities that make perturbation theory inconsistent. 

The purpose of the first part of this paper (Sec. 2) is to consider the condition 
for the perturbative QCD coupling to be causal as a function of Nf (for a general 
N c ), and compare it with the lower boundary of the conformal window set by su- 
perconvergence. In Sec. 2.1 we study causality at the level of the 2-loop coupling 
and in Sec. 2.2 we examine the effect of higher orders. Next, in Sec. 3, we study 
the same issue in SQCD. In Sec. 3.1 we investigate the singularity structure of the 
2-loop coupling in the electric theory and study the effect of higher orders. In Sec. 
3.2 we discuss the dual (magnetic) theory. Sec. 3.3 summarizes the main findings of 
Sec. 2 and 3. In Sec. 4 we consider the Banks-Zaks expansion for the value of the 
fixed-point (Sec. 4.1) and for the critical exponent 7 that controls the rate at which 
a generic effective charge approaches the fixed-point in the infrared limit (Sec. 4.2). 
In Sec. 4.3 we show that the 2-point Pade approximants technique [21] can be used 
to interpolate between the Banks-Zaks expansions for a physical quantity in the two 
dual theories. The example considered is the critical exponent in the large N c limit. 
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2 The conformal window in QCD and the analyt- 
icity structure of the coupling 

As explained in the introduction, according to the superconvergence criterion [4, 5] 
an SU(N C ) gauge theory with Nf light flavors is confining so long as the anoma- 
lous dimension of the gauge field propagator 7 00 of eq. (4) is negative, i.e. for 
r = Nf/N c < 13/4. This is only a sufficient condition for confinement, and there- 
fore we can expect a phase transition from the confining phase to a phase which is 
conformally invariant in the infrared, either at r = 13/4 = 3.25 or somewhere above 
this line. 

Referring to the superconvergence criterion as determining the lower boundary 
of the conformal window, we now turn to study the perturbative f3 function. In 
real-world QCD the perturbative running coupling has "causality violating" Landau 
singularities on the space-like axis which, according to the common lore, signal the 
inapplicability of perturbation theory in the infrared region and the necessity of non- 
perturbative power like terms. On the other hand, close enough to the top of the 
conformal window causality can be established within the perturbative framework 
[11, 12]. There the perturbative (3 function leads to a causal running coupling, which 
has a finite infrared limit and no Landau singularities in the entire Q 2 plane: its 
only discontinuity is a cut along the time-like axis. In this situation the perturbative 
analysis does not signal the need for non-perturbative physics. It is then possible 
that perturbation theory by itself describes well the infrared physics. 

By definition, in the conformal window the coupling reaches a finite limit in 
the infrared. As explained above, in the upper part of the window this finite limit 
is obtainable from the perturbative f3 function. Is it true also away from the top 
of the window? In other words, is it the perturbative coupling that reaches a finite 
limit? and in this case, can we reliably calculate the fixed-point value in perturbation 
theory? In order to address these questions we study here the conditions for a causal 
perturbative coupling and compare them to the boundary of the conformal window. 

2.1 Causality from the 2-loop (3 function 

The 2-loop (3 function with (3q > and (3\ < is the simplest example where Landau 
singularities can be avoided, so it is natural to begin by analyzing this example. It 
should be stressed that the 2-loop (3 function corresponds to a particular choice of 
renormalization scheme, the so-called 't Hooft scheme, where all the higher-order 
corrections to the (3 function (3 2 , 03 and on vanish. Since the first two coefficients of 
the f3 function are scheme invariant, we shall obtain a criterion for causality which 
does not have an explicit dependence on the scheme. 
The 2-loop renormalization group equation, 




(6) 
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where t = ln(Q 2 /A 2 ), and [3] 



c = 



Pi J_ 
Po 16^o 



34 *r2 



i 



A 7 



(7) 



can be integrated exactly [12] using the Lambert W function [22]. It was shown 
in [12] that if c > a Landau branch point is present on the space-like axis, if 
— Po < c < a pair of Landau branch points appears at some complex Q 2 values 
and if 

c < -p < (8) 

the coupling has a causal analyticity structure, with no Landau singularities^. We 
rederive here these results using a simpler (but less rigorous) approach [23]. 
Integrating (6) we obtain 



ln(Q 2 /A 2 



1 



P X ^-loop 



In 



1 

X-p-p 



(9) 



where x PP = — 1/c and 7 2 loop is the critical exponent at the 2-loop order. 7 is 
defined as the derivative of the ft function at the fixed point, 



dp(x) 

7 = 



dx 



(10) 



x=xpp 



At 2-loop order 7 2 -'°°p = -p 2 jp v . 

If c > there is a Landau singularity on the space-like axis. A positive fixed- 
point is obtained for c < 0, but this condition alone does not guarantee causality - 
there can be Landau singularities in the complex Q 2 plane. 

Assuming that the singularities are such that |x(<5 2 )| — > 00, we expand (9) 
around these points in powers of 1/x. The leading term in this expansion gives the 
location of the singularity. The phase of the r.h.s. of (9) at the singularity is 



±7T 



1 



7 



2-loop ' 



(11) 



If |$| < 7r, i.e. if 7 2 loop > 1, or — p < c < 0, the singularities are in the first sheet 
(the time-like axis cut corresponds to |<&| = vr), and if |$| > 7r, i.e. if 



< 7 



2— loop 



< 1 



;i2) 



or c < — Po, the 2-loop coupling is causal. The condition (12) (or (8)) for a causal 2- 
loop coupling translates in QCD into the following condition for r upon substituting 
Po and c from eqs. (2) and (7), respectively: 



4r 2 + 



9 



r + 223 < 



(13) 



§ We always assume asymptotic freedom, i.e. /3 > 0. 
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This leads to an approximately N c independent critical value for r for any possible 
value of N c (since iV c 2 3> 9/83), namely the 2-loop coupling is causal as long as 

r> (83-9-n/41)/8~3.17. (14) 

For lower r, the condition (12) does not hold and there appears a pair of complex 
singularities in the Q 2 plane. If r is reduced further, c becomes positive and then a 
Landau branch point appears on the space-like axis. This change occurs at: 

~ 2.62. (15) 

The results are summarized in fig. 1 in the upper plot, where the lower boundary 
of the conformal window implied by superconvergence (r = 3.25) is compared with 
the lower boundary of the region where the 2-loop coupling is causal according to 
(13), which is asymptotic at large N c to r ~ 3.17. Clearly, the 2-loop coupling is 
causal in the entire conformal window. This conclusion holds, of course, also in the 
case where the lower boundary of the conformal window is somewhere above the 
critical value for superconvergence (r = 3.25). This suggests that the fixed-point in 
QCD is always of perturbative origin. 

The proximity of the two lines, the upper boundary of the superconvergence 
region (r = 3.25) and the lower boundary of the 2-loop causality region (r = 3.17) 
does not have any deep meaning. Presence of complex Landau singularities in the 
running coupling signals that the coupling becomes strong but it does not necessarily 
imply confinement - an example is provided by SQCD (Sec. 3). 

Due to the closeness of the two lines one might worry that even within the confor- 
mal window the large distance physics cannot be reliably described by perturbation 
theory. However, we shall see in the next section that 3-loop corrections make the 
coupling causal in a wider range, and eventually perturbation theory does seem 
reliable down to the bottom of the conformal window. 

2.2 How relevant is the criterion for causality at 2-loop? 

It is natural to wonder whether the singularity structure of the coupling which is 
defined by the truncated 2-loop f3 function is of any physical significance. Of course 
we do not doubt the assumption that the theory as a whole is causal. According 
to the common lore, the appearance of the non-physical Landau singularity in the 
perturbative coupling in real-world QCD is nothing more than a sign of the inappli- 
cability of perturbation theory for describing the infrared region. Thus the presence 
of Landau singularity indicates the significance of non-perturbative corrections in 
the infrared. The interesting point is that close enough to the top of the conformal 
window, there may be a possibility to establish causality using only perturbation 
theory, as we explain below. 



r = 34/ 
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2.2.1 Causality beyond 2-loop - general discussion 



In general, the analyticity structure of a coupling based on some higher order (3 
function, 

P(x) = -p x 2 [l + cx + c 2 x 2 + ■ ■ •] (16) 

can be completely different from that of the 2-loop coupling (6). This is clearly so if 
Landau singularities are present: their location and nature generally depend on all 
the coefficients of the /3 function and consequently on the renormalization scheme. 
This "instability" should be of no surprise since the weak coupling expansion breaks 
down completely when examining the singularities of the coupling. 

As an example how the singularity structure changes and becomes more complex 
as higher order terms in the (3 function are included, consider the 1-loop coupling, 
the 2-loop coupling and Pade improved 3-loop coupling, defined by 

R <*.\ a [c-(c 2 /c)]x 

which were all analyzed in [12]. The 1-loop coupling has a space-like Landau pole, 
the 2-loop coupling can have a causal structure or a pair of complex branch points 
or a space- like branch point. The Pade improved 3-loop coupling can be causal but 
it can also have both simple poles and branch points (the details appear in [12]). 

While these examples show that there is no stability when going to higher orders 
if Landau singularities exist, they also indicate that if the 2-loop coupling is causal, 
causality may be preserved when higher order corrections are included. In fact, it 
is rather simple to explain why this kind of stability be expected in general. When 
the 2-loop coupling is causal it is bounded, and in many cases also small, for any 
complex Q 2 . If so the usual perturbative justification holds: the next term in the (3 
function series which is proportional to a higher power of the coupling is small, and 
likewise higher order terms. In this situation higher-order terms are not expected to 
have much influence on x{Q 2 ). In other words, absence of Landau singularities can 
be consistently confirmed at the perturbative level, whereas presence of Landau sin- 
gularities can only be confirmed or disproved in the full theory by non-perturbative 
methods. 

The first step in establishing causality in perturbation theory is to examine the 
analyticity structure of the 2-loop coupling, as we did above. On one hand the 2-loop 
coupling has the advantage that it does not depend on the renormalization scheme. 
On the other hand it does not correspond directly to any observable and therefore 
it may not be causal. Thus, we are forced to examine higher order corrections 
(or renormalization schemes other than the 't Hooft scheme), and see whether the 
causality condition at 2-loop order is reasonable. The next step is therefore to 
choose a representative renormalization scheme, different from the 't Hooft scheme, 
and ask whether the 3-loop correction to the (3 function has a significant effect on 
the infrared coupling. If the effect of the 3-loop correction is negligible, that is if 
the coupling is small enough such that 

(3 2 x 2 {Q 2 ) « (3 lX (Q 2 ) (18) 
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in the entire complex Q 2 plane, then we shall consider that causality is established 
at the perturbative level. 

To be completely convinced, one might want to check also the magnitude of 
higher-order corrections corresponding to 4-loop order and beyond. However, it is 
important to remember in this respect, that if we go to high enough order (n), we 
will always obtain 

\P n x n \ > \f3 n . lX n - l \ (19) 

due to the asymptotic nature of the (3 function series, and thus it does not make 
sense to require that all the higher-order terms will be small. In the scenario de- 
scribed above, namely that the 2-loop coupling is already causal and small it seems 
reasonable to require that the 3-loop correction is small and stop there. Clearly, 
this scenario is just the simplest case to consider. It is possible that x(Q 2 ) at 2-loop 
is still not small enough so as to guarantee \(3 2 x 2 \ <C \f3\x\, but f3 2 is negative so 
x(Q 2 ) at 3-loop is much smaller, and then higher order corrections are negligible: 
\03X 3 \ <C |/92^ 2 |- Of course, in this case the results might depend on the renormal- 
ization scheme. 

An encouraging observation with regards to the 2-loop analysis is that the con- 
dition for causality of the 3-loop coupling is quite modest once the 2-loop coupling is 
causal. We will show that the only further requirement is that the 3-loop (5 function 
has a positive root corresponding to the infrared stable fixed-point. 

It is most convenient for this demonstration to write the 3-loop f3 function in the 
following form: 

m = i^m = I /i/2 _ (/i + h)x + x2 l (20) 

where fj 2 = (3 /(3 2 and f 1 + f 2 = — /?i//3 2 , thus: 

A = ^(-ft + VA) (21) 
h = ^(-ft-VA) 



with A = j3f — 4:f3 f3 2 . Note that in all cases of interest, namely when there is a 
positive real zero to the 3-loop (5 function, the infrared fixed-point is x m = f 2 > 0, 
and x vv = fi is an ultraviolet fixed-point. The corresponding critical exponents are, 
for x m : 

l!n l °° P = P2X 2 m (x vv - x m ) = x 2 m VX > (22) 

and for x vv : 

7uv °° P = ~P2X 2 vv {x vv - x m ) = -x 2 vv VA < 0. (23) 



It is useful to note that 



+ 1 = 1 = (24) 

~ 3-loon -9,-lnn-n Oil \ > 



3-loop 3-loop ^2-loop ' /32' 

luv Jm I 
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where we have used the definition of f\ : 2- It then follows, assuming 2-loop causality 
(12), that 

< 7?r l °° p < l 2 ~ loop < 1. (25) 
In order to examine causality at 3-loop we integrate (20) and obtain: 



HQ 2 /A 2 ) = ^ + ^ln 

Pox 7 UV 



1 



Ob 0b\ 



1 . 

7m 



1 1 

X x m 



(26) 



To find the causality condition, we study, as in the 2-loop case, the phase of the 
Landau singularity 1 . We assume that the only singularities are such that \x{Q 2 ) | — > 
oo, and expand (26) around these points in powers of 1/x. The leading term in this 
expansion gives the location of the singularity If fa < 0, then x uv is negative and 
the phase of the r.h.s. of (26) at the singularity is 

7ir 

while if > 0, x VY is positive and the phase is 

( 3— loop 3-loop ) 71 ^2-loop ' (^S) 

7uv 7ir / " 

Using (25) we find that in both cases |<E»| > ir and it follows that the 3-loop coupling 
is causal. 

We showed that if the 2-loop coupling is causal, the 3-loop coupling is also causal, 
provided it has an infrared fixed-point. It is not clear whether such a conclusion 
can be extended to higher orders. It may be however interesting to note that we 
already know from the analysis of [12] another example where similar conclusions 
hold: this is the Pade improved 3-loop j3 function, defined by (17). Contrary to 
the above examples, this (3 function is not truncated at some finite order, and thus 
it could be expected a priori to behave differently. According to [12] the causality 
condition for the Pade- improved 3-loop coupling is c < — (3 and c 2 < c 2 . The first 
condition is the same as the condition for the causality of the 2-loop coupling. In 
fact, the critical exponent of this coupling is equal to that of the 2-loop coupling, 
and thus the condition is < 7 PA = '-f 2 ~ loop < 1. The second condition is just the 
condition to have a positive infrared fixed-point. 

We comment that the inverse statement does not hold: the 3-loop coupling can 
be causal even if the 2-loop coupling is not. If (3 2 is negative and large enough, the 
3-loop coupling is causal independently of the sign of f3\. 

Coming to analyze the conditions for causality or the stability of the causal 
solution with respect to higher order corrections (such as (18)), we should, in general 
solve the renormalization group equation at each order to obtain x{Q 2 ) in the entire 
Q 2 plane, as was done in [12] for the 2-loop and the Pade improved 3-loop couplings. 

1 As opposed to the 2-loop case, where it is also possible to invert [12] the relation (9) to calculate 
x(Q 2 ) using the Lambert W function, here this cannot be done. 
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However, we shall demonstrate below that it is in fact enough to examine the effect 
of higher orders on the infrared limit of the space-like coupling x PP = x(0), unless 
the coefficients of the f3 function are extremely close to the condition where causality 
is lost. In most cases when x{Q 2 ) is causal, |rc((5 2 )| ^ x(0) in the entire Q 2 plane. 
Of course, when causality is lost x{Q 2 ) diverges at some point while x(0) is finite, 
and thus close to the boundary of the causality region x(0) is not indicative at all. 
The point is that the region where the maximum value of |x(Q 2 )| is much larger 
than x(0) is quite narrow. To demonstrate this, consider again the example of 2- 
loop coupling in N c = 3 QCD. According to (13) this coupling is causal so long as 
Nf k, 9.683. At the point where causality is lost, \x(Q 2 )\ reaches infinity on the 
first sheet (on the time-like axis), while the space-like coupling has its maximum at 
x(0) ~ 0.88 which is not so large. However, if we move slightly above the causality 
boundary, the maximal value of |x((5 2 )| in the entire Q 2 plane becomes of the order 
of x(0). We show this phenomenon in fig. 2 where we plot the region in the complex 
coupling plane, into which the entire complex Q 2 plane is mapped. The contour 
itself corresponds to the cut along the time-like axis (Q 2 < 0) and it was computed 
using the Lambert W function solution, as explained in [12, 11]. As shown in the 
plot, for Nf = 9.7, i.e. very close to the point where causality is lost, the maximal 
value of |x(<5 2 )| on the time-like axis is still significantly larger than x(0). One 
clearly identifies here the effects of the singularities that are present on the second 
sheet. On the other hand, already at Nf = 10, the maximal value of |a;(<3 2 )| on the 
time-like axis is of the order of x(0). 

We found that the condition for the 2-loop coupling to have a causal analyticity 
structure is < ^ 2 ~ loop < 1, where ^ 2 ~'°°p = corresponds to a free theory, the 
limit obtained at the top of the conformal window, and ^ 2 ^ 1o °p = 1 corresponds to 
the point where Landau singularities first appear. At 2-loop order the condition 
< '-f 2 ~ loop < 1 is both a sufficient and a necessary condition. It is interesting to 
see how this generalizes to higher orders. When the j3 function has more than one 
zero, we should specify at which of them 7 is defined. The only root that is relevant 
in the asymptotically free phase is the smallest positive zero, the physical infrared 
stable fixed-point, and we always refer to this one. 

The 3-loop analysis shows that < ^ 3 ~ 1o °p < 1 is a necessary condition but not 
a sufficient one. An example where the 3-loop coupling is not causal although the 
above condition is obeyed can be constructed starting with a non-causal 2-loop (5 
function with fa < and adding a 3-loop term with fa positive but small enough 
such that a positive zero for the 3-loop (5 function exists. It then follows from eq. 
(28) that the 3-loop coupling is not causal although ^ 3 ~ 1o °p can still obey the above 
condition. We stress that this example is not representative since usually, as we 
shall see, fa < and then the condition < 7 < 1 is both necessary and sufficient 
also at the 3-loop order. 

In fact the condition < 7 < 1 is always necessary for a causal analyticity 
structure. The condition 7 > is simply the one to have an infrared stable fixed- 
point. To show that also 7 < 1 is necessary we use the following observation: a 
causal structure implies that there is a well defined mapping x{Q 2 ) from the entire 
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complex Q 2 (the first sheet) into a compact domain in the complex coupling plane, 
such that for large enough \Q 2 \ the coupling flows to the trivial fixed-point, as 
implied by asymptotic freedom". As we saw in the example of fig. 2 (these features 
as completely general) the space-like axis is mapped to real positive values in the 
range [0, x FP ] and the time-like axis is mapped to the boundary of this domain in the 
complex coupling plane. It follows from the definition of 7 in (10) that the coupling 
approaches the fixed-point according to 

i=irp ~(S) 7 ' (29) 

where A cff is an observable-dependent QCD scale. If 7 > 1, there is a phase $ = ir/'-f 
in the complex Q 2 plane (Q 2 = exp(i<I>)) such that in the limit Ql — ► the rays 
±$ are mapped by (29) to positive real values of the coupling larger than the 
fixed-point value (x = x+ p ). On the other hand a straightforward analysis of the 
(3 function shows that values of the coupling x > x FP either belong to the domain 
of attraction of some non-trivial ultraviolet fixed-point or flow to an ultraviolet 
Landau singularity. The conclusion is that there is no singularity free mapping that 
obeys the asymptotic freedom condition stated above. In particular, if two different 
ultraviolet fixed-points are allowed for different values of Q 2 it implies the existence 
of a separatrix, discriminating between the values of Q 2 that flow to each of the 
ultraviolet fixed-points, i.e. there are singularities in the first sheet of the complex 
Q 2 plane. We stress that the arguments why < 7 < 1 is a necessary condition 
for a causal analyticity structure are completely general: they are not based on 
perturbation theory. 

2.2.2 Causality at higher orders in QCD 

We would like to examine whether causality can be established in perturbation 
theory in the specific case of the conformal window in QCD. Close to the top of the 
conformal window, causality is established at the 2-loop level. The infrared coupling 
is small and thus the 3-loop term is negligible and condition (18) for stability of the 
perturbative analysis is obeyed. This is no longer true at the bottom of the window. 

We start the discussion in the MS scheme, which has the advantage that the 4- 
loop coefficient in the (3 function is known [24, 25] . We shall refer to physical effective 
charges later. The fixed-point value of the coupling, calculated as an explicit solution 
of the equation /3(x) = in the large N c limit at 2-loop and then in the MS scheme 
at 3-loop and 4-loop orders is shown in fig. 3 as a function of the distance from the 
top of the conformal window, 

11 11 N f 

-r = -f. (30) 



2 2 N, 



c 



Hit was demonstrated in [12] in the particular case of the Lambert W solution for the 2-loop 
coupling that in order to define the analytical continuation of x(Q 2 ) from the space- like axis to 
the entire first sheet, it is essential to require asymptotic freedom for complex Q 2 values. 
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The 2-loop coupling reaches relatively large values towards the bottom of the con- 
formal window, but then the 3-loop and 4-loop couplings take significantly lower 
values, and in addition they are very close to each other. These results can be un- 
derstood knowing the negative sign of the 3-loop coefficient fa and the magnitude 
of successive terms in the MS scheme, shown in fig. 4. In the latter, the coupling 
is evaluated at the fixed-point according to the zero of the 3-loop (5 function. It is 
clear from the plot that the condition for stability of the 2-loop result (18) does not 
hold in the lower part of the conformal window. It certainly does not hold for r ^ 4, 
corresponding to Nf ^ 12 since there the 3-loop term is comparable to the 2-loop 
term. Thus we are forced to examine causality at higher orders. 

Since the 3-loop coefficient in MS is negative for the relevant Nf/N c values, 
the 3-loop (5 function has a positive real fixed-point, and according to the general 
discussion in the previous section, 3-loop causality is implied within the region where 
the 2-loop coupling is causal. Now, in order to trust the 3-loop causality, it is 
required that the 4-loop term will be small enough. Indeed, as shown in fig. 4 the 
4-loop term in the MS scheme remains small in the entire conformal window. The 
effect of the 4-loop term on the fixed-point value is shown in fig. 3. Clearly, this 
is a negligible effect, and thus perturbative stability is realized at the 3-loop level. 
It would be better to check the effect of the 4-loop term on x(Q 2 ) in the entire 
Q 2 plane, but based on the experience with the 2-loop coupling we expect that in 
general the space-like fixed-point value is indicative of the magnitude of x{Q 2 ) in 
the entire complex Q 2 plane. 

Next we consider the value of the critical exponent as a function of the distance 
from the top of the window. The results of an explicit calculation of 7, in the 
large N c limit, from the 2-loop, 3-loop and 4-loop f3 functions in MS are shown 
in fig. 5 in the upper plot. In agreement with our previous discussion the condi- 
tion < 7 3_/oo p < ^-loop < 1 j s obeyed in the entire conformal window. The points 
where the 2-loop and 3-loop couplings cease to be causal can be identified in this 
figure as the points where 7 = 1. Since the 4-loop term is small, ^ 4 ~ loop ~ ^~ 1oo p 
within the resolution of this plot, and so the perturbative stability which charac- 
terises the coupling exists also for the critical exponent. 

We stress that the results described above are not special to the large iV c limit. 
In particular fig. 3 through 5 are qualitatively the same for any iV c . 

The above investigation shows that the MS coupling is causal at the 3-loop level 
and, given the smallness of the 4-loop term, presumably also at the 4-loop level in 
the entire conformal window. However, this coupling does not correspond directly 
to any observable quantity. It is important to check whether similar conclusions 
apply in physical schemes. 

A relevant analysis has been performed in [11]. Fig. 1 in [11] compares the Nf 
dependence of c 2 = fa/ 'fa for various physical effective charges. The observation 
that c 2 for different effective charges are numerically close and that they share the 
same Nf dependence indicates that certain properties of the coupling may be generic 
in spite of scheme dependence. In particular we note that c 2 is negative in the entire 
conformal window not only in MS , but also for all the physical effective charges 
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considered. We conclude that there is a fixed-point at the 3-loop order in all these 
physical schemes and, according to the general discussion above, 3-loop causality 
follows**. Unfortunately, 4- loop coefficients in physical renormalization schemes are 
not known yet^. Consequently, the stability of the 3-loop causal coupling with 
respect to higher loop corrections cannot be studied for physical effective charges 
like we did in the MS scheme. However an alternative is provided by the Banks-Zaks 
expansion, which can be calculated in physical schemes up to next-to-next-to-leading 
order term [14, 15, 11]. This will be discussed further in Sec. 4. 

We comment that the perturbative coupling at the 3-loop order can have a causal 
analyticity structure even somewhat below the bottom of the conformal window, i.e. 
in the upper part of the confining phase. In this respect, different couplings may 
behave differently. We recall that for N c = 3 the 't Hooft coupling, defined by the 
truncated 2-loop f3 function, is causal down to Nf ~ 9.68, quite close to the bottom 
of the conformal window Nf ~ 9.75. This can be compared with the MS scheme 
where 3-loop causality is lost at Nf ~ 8.5 and to physical renormalization schemes 
in which the causality domain is even wider. Based on the results of [11] and the 
above type of analysis we find that for the effective charge defined by the vacuum 
polarization D-function and the ones associated with the polarized and non-polarized 
Bjorken sum-rules 3-loop causality is lost at Nf ~ 7.2 while for the effective charge 
defined from the heavy quark potential 3-loop causality is lost at Nf ~ 8.4 H . 

Finally we consider the calculation of the critical exponent using physical renor- 
malization schemes. As long as the fixed point is perturbative, it is natural to expect 
that 7 could be calculated with a reasonable accuracy starting with the truncated 
f3 function in various renormalization schemes. Since 7 is a universal quantity the 
results should agree. The results of an explicit calculation of 7 in several physi- 
cal schemes at the 3-loop order are presented in fig. 5 together with the results in 
MS . The schemes we use include the vacuum polarization D-function, the polarized 
and non-polarized Bjorken sum-rules (the latter two curves overlap) and the heavy 
quark effective potential. The results in the different schemes agree very well close 
to the top of the window. The spread increases to about ±15% towards the bottom 
of the window and is interpreted as an artifact of using a truncated perturbative 
expansion. We shall come back to discuss the accuracy to which 7 can be calculated 
in sec. 4.2.1 in the framework of the Banks-Zaks expansion (see table 5 there). 

**Note that the Pade improved 3-loop coupling is also causal. 

tt An exception is the effective charge related to the Higgs hadronic decay width. For this quantity 
the infrared fixed-point does not even exists in the lower part of the conformal window due to a 
large positive 4-loop coefficient. We do not, however, consider this example as representative (see 
the discussion in [11]). 

ttThe last result is based on the recently published 2-loop calculation of the static potential in 
QCD [26], which corrects a previous result used in [11]. 
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3 The conformal window in SQCD 



The (3 function in SQCD is given by* 

P(X) = ~T7~~T7y2\ = "(^ 2 + B lX 3 + •••) = -B^X 2 (1 + C lX + ■ ■ •) (31) 

where x = a /it = g 2 /(An 2 ), 

B = \ (3iV c - N f ) , (32) 

and 

Cl S * = l -N c !Lt!, (33) 

1 B 2 c 3N C -N f 2N C ' 1 ; 

where the coefficients where calculated in [28]. Above the line i? = Nf/N c = 3 the 
theory is infrared free, while below this line it is asymptotically free in the ultraviolet. 
For R just below 3, B is small and positive and B\ is negative, leading to an infrared 
fixed-point at a small x FP ~ —B Q /Bi = — 1/Ci, making the theory weakly coupled 
at all scales. As Nf (and thus R) is decreased, the infrared coupling increases. 
According to Seiberg [16], the infrared fixed-point persists even down to such low 
Nf that the original degrees of freedom are strongly coupled and then a dual theory 
which is based on another gauge group with Nf = Nf — N c colors is appropriate to 
describe the infrared limit (d stands for a dual variable). Seiberg's conjecture can 
only be understood if the fixed-point is of non-perturbative origin, at least in the 
lower part of the conformal window. This is contrary to our previous observation 
concerning the perturbative origin of the fixed-point in the non-supersymmetric case. 
Thus we would like to check that indeed a definite difference exists between the 
conformal window in QCD, which is perturbative and the one in SQCD which is 
not. This is done here by considering the analyticity structure of the coupling 
constant and in the next section, by comparing the Banks-Zaks expansion in SQCD 
to that in QCD. We shall indeed see that already at the perturbative level SQCD is 
more strongly coupled than QCD in the lower part of the conformal window. The 
fact that the strong coupling nature of SQCD at the lower part of the window is 
manifested in perturbation theory is not obvious a priori; strong infrared effects 
could have been induced instead by terms invisible to perturbation theory. 

Duality [16] provides an intuitive description of the conformal window in SQCD, 
which is absent in QCD. The lower boundary of the conformal window in SQCD 
is naturally identified as the R = Nf/N c ratio at which the dual theory undergoes 
a phase transition from the asymptotically free phase (inside the window) to the 
infrared free phase (below the window). The 1-loop (5 function coefficient in the 
dual theory can be obtained by substituting Nf — N c for N c in (32): 

B d = \ (2N f - 3N C ) . (34) 



* Capital letters are used here to distinguish SQCD coefficients from QCD ones. 



15 



Thus, the conformal window is3/2<i?<3, as shown in the lower plot of fig. 1. 
The original theory is weakly coupled, and therefore provides a natural physical 
description, for R just below the line R = 3, while the dual theory is weakly coupled 
just above the line R = 3/2. 

An important consistency check for both duality and the superconvergence cri- 
terion for confinement is that the lower boundary of the conformal window in both 
approaches coincides [19]*. The observation of [19] is the following: in SQCD the 
anomalous dimension of the gluon propagator in the Landau gauge is 

7o S o QCD = ~\ (|* c - N f ) , (35) 

which is just proportional to the first coefficient of the /3 function Bq in the dual 
theory. As a result, 7o(? CD becomes negative, implying superconvergence and there- 
fore confinement for the original theory, as R becomes smaller than 3/2, i.e. exactly 
where the dual theory becomes infrared free (Bq in (34) changes sign). 



3.1 The analyticity structure of the SQCD coupling 

The purpose of this section is to analyze the singularity structure of the perturbative 
SQCD coupling, in parallel with the analysis of the QCD coupling in Sec. 2, and in 
particular to find when it is consistent with causality. 

The first step is to analyse the 2-loop coupling. The 2-loop causality condition 
C\ < —Bq translates, using eqs. (32), and (33), to the following condition for R: 



R 2 + 



2 

m~ 10 



i? + 15<0. (36) 



Similarly to the non-supersymmetric case (see eq. (13)), the condition (36) leads to 
an approximately N c independent critical value for R for any possible value of N c 
(since N% 3> 2/10), namely the 2-loop coupling is causal as long as 

R> 5-^- 1.8377. (37) 

The crucial observation is that the line (37) that limits from below the region 
where the 2-loop coupling is causal, is within the conformal window which has its 
lower boundary at R = 3/2. This is shown in the lower plot of fig. 1. The situation 
encountered here is contrary to the one in non-supersymmetric QCD, where the 
2-loop perturbative coupling is causal in the entire conformal window. This obser- 
vation fits the general expectation based on duality, that the fixed-point in SQCD 
is non-perturbative in the lower part of the conformal window. 

In addition we ask when does the 2-loop coupling develop a space-like Landau 
singularity. The condition C 1 > translates (using (33)) to the following: 

* < ^wwr (38) 

'''The generalization of this result to other supersymmetric models was examined in [20]. 
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which is asymptotic in the large N c limit to the lower boundary of the conformal 
window, R = 3/2. 

Note that in the supersymmetric case, it is natural to use the NSVZ form [28, 29] 
of the (3 function, and thus one may wonder if our results concerning the analyticity 
structure of the coupling may vary when using the (3 function in this form rather than 
the truncated 2-loop one. In the Appendix we show that both the condition for a 
causal coupling and the condition for a space-like Landau singularity are exactly the 
same in the two cases, if in the NSVZ form one uses the leading order approximation 
for the matter field anomalous dimension. 

The next step in the analysis of the perturbative coupling causality, as in the 
QCD case, should be to examine the effect of higher order terms in the perturbative 
(3 function. We choose to work in the DRED renormalization scheme [27], assuming 
that our conclusions will not depend on this choice. 

The explicit solutions of (3{x) = for A^ c x(0) in the large N c limit are shown in 
fig. 6 as a function of the distance from the top of the conformal window, 

5 = 3- R = 3- -L. (39) 

The 2-loop solution is infinite at the bottom of the window (see (38)). Already here 
we encounter a situation different from QCD, namely stronger coupling. Since the 
3-loop coefficient is negative, the 3-loop solution is smaller. The latter is finite down 
to the bottom of the window, but it is still rather large. The fixed-point at 4-loop 
order exists only up to 5 ~ 0.4 (near the 4-loop arrow in the figure). Beyond this 
point there is no positive real solution to the equation (3{x) = 0. The reason is that 
the 4-loop term in SQCD is positive (like in QCD) and large (contrary to QCD) 
as can be learned from fig. 7. This figure shows the relative magnitude of the four 
leading terms in the large N c SQCD (3 function. The coupling in fig. 7 is evaluated 
as the zero of the 3-loop f3 function. 

In the lower plot of fig. 5 we show the value of the critical exponent as a function 
of 8 according to the 2-loop, 3-loop and 4-loop order large N c DRED f3 function. 
The necessary condition for a causal structure 7 = 1 is reached by both the 2-loop 
order, which was discussed above, and 3-loop order solutions for 7 well within the 
conformal window. The 4-loop result for 7 exists of course only up to 5 ~ 0.4 where 
a positive fixed-point exists. 

Examining fig. 5 through 7 we can determine where causality can be established 
in SQCD at the perturbative level. In the upper part of the conformal window the 
3-loop term is small with respect to the 2-loop one, so one can trust 2-loop causality. 
As R is decreased the 3-loop term becomes comparable to the 2-loop term and then 
one has to consider causality at 3-loop order. The negative sign of fa guarantees 
that the 3-loop coupling is causal at least as long as the 2-loop coupling is. But since 
the 4-loop term is very large, the perturbative argumentation fails. Thus in SQCD, 
it is possible to establish causality in perturbation theory only in the upper part of 
the conformal window. To be specific, two alternative criterions can be considered: 
the first is to require that the 3-loop term will be smaller than the 2-loop term. The 
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two become equal around R ~ 1.9, i.e. just above the 2-loop causality boundary. 
The second is even more restrictive, namely to require that also the 4-loop term is 
small, or that the 4-loop (3 function will have a positive real root. This is realized 
only above R ~ 2.6. 

Maybe the most interesting observation in fig. 7 is the fact that the 4-loop term 
in the SQCD (3 function is larger than the leading terms already very close to the 
top of the conformal window. This may be related to the asymptotic nature of the 
(3 function series. The asymptotic behavior is another aspect in which the SQCD f3 
function is presumably different from the QCD one, a point which certainly deserves 
further study. 



3.2 Reduction of Couplings in the magnetic theory 

In the previous section we studied the singularity structure of the 2-loop coupling 
in the electric theory. Our aim here is to perform a parallel analysis in its dual, the 
magnetic theory. This is, however, not straightforward since the magnetic theory 
has two couplings, rather than one. The running of the gauge coupling is affected 
by the Yukawa interaction of the chiral quark superfields with the mesons, which 
is described by the superpotential (5). This gives rise to coupled renormalization 
group equations of the form 

£{ X ,\) = -J?— = -B d x*-B d x*-B d x*\ + --- (40) 



cHn(Q 2 ) 



with 



B d = \(2N f -3N c ) (41) 



8 v 1 1 8 N f -N c 

1 

I' 



C 



d _ HN f -N c ) 2 



A 



2 N f -N c 
C d xx = ±(3N f -N e ) 

where Bq and Bf can be obtained by substituting N c — > Nf — N c in Bq and Bi, and 
the other coefficients where calculated in [30]. Note that in (41) we use N c to denote 
the number of colors in the original (electric) theory, and thus the dual theory has 
an SU (Nf — iV c ) gauge symmetry. This is contrary to the notation used in [30] that 
corresponds to an SU(N C ) gauge group in the magnetic theory. In addition, note 
that in [30] there is a typo in eq. (64), where a factor of two is missing in the second 
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term in the second equation*. The correct factor can be easily obtained by using 
eqs. (55), (62) and (63) there. Our coefficients do agree with those in [31]. 

In order to study the analyticity structure of the coupling in the magnetic theory 
one should, in principle, integrate the coupled renormalization-group equation (40). 
This is, however, rather complicated, and so we choose a simpler approach (which 
remains to be further justified) based on the notion of Reduction of Couplings. 

It was recently shown by Oehme [31] that there is a unique reduction of the 
coupled renormalization-group equation (40) to a single-coupling equation such that 
the superpotential does not vanish, which is essential for duality. Ref. [31] describes 
in detail how to apply the general method of Reduction of Couplings to this problem. 
We shall use here only the leading order relation between A and x. To obtain the 
relation between the couplings one assumes 

\(x)=f(N c ,N f )x + 0(x 2 ) (42) 

and imposes the consistency condition, 

(3 d (x,\(x)) = ^(3 d (x,\(x)). (43) 
Using (40), the condition (43) leads, at leading order, to: 

f(N c , N f ) [C d xx f(N c , N f ) + C d + B d ] = (44) 
and for a non-vanishing superpotential, the results is 

/liVc ' iV/j C d xx 2(N f -N c )(3N f -N c y ^ } 

Note that f(N c ,Nf) is positive in the entire conformal window. 

With the result (45) at hand we can substitute the A term for f(N c , Nf) x in the 
equation of (3 d (x, A) and obtain a single-coupling renormalization-group equation 
which is valid up to 2-loop order: 

^ x{x)) = = ~ Bd ° x2 - B i x3 +--- ( 46 ) 

where 

B d — Bf + Bf x f(N c , Nf). (47) 

Next, we would like to analyze the analyticity structure of the coupling in the 
dual theory, using the reduced f3 function (46). Let us calculate first the condition 
for the dual 2-loop coupling to have a causal analyticity structure (the analog of 
(37) in the original theory). The causality condition, Bf/B^ < —Bq, yields 



18i? 3 - 6AR 2 + 



65 + 4 



R - 15 < (48) 



*The authors thank D. Anselmi and R. Ochmc for their help on this matter. 
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which again leads to an approximately N c independent critical value for R for any 
possible value of N c (since N 2 3> 2/65), namely the 2-loop coupling in the dual 
theory is causal as long as 

R & 1.8357. (49) 

As with the original theory, the 2-loop causality region of the dual perturbative 
coupling does not cover the far-end of the window. Note (fig. 1) that the regions of 
a causal 2-loop coupling in the two dual descriptions, (37) and (49) do not overlap. 
This fits the intuition on which duality is based, i.e. that when one theory is weakly 
coupled its dual is necessarily strongly coupled. Since we assume that within the 
window a consistent perturbation theory implies small non-perturbative effects, an 
overlap would lead to contradiction: it would suggest that two different weakly 
coupled theories can describe the same infrared physics. 

In fig. 1 the 2-loop causality boundaries in the two theories are very close. How- 
ever, if one adopts a conservative attitude that perturbation theory actually breaks 
down above the 2-loop causality boundary (taking into account the large 4-loop 
correction) the perturbative regions of the two theories will be more separated. 

One can also find the condition to have no space-like Landau singularity in the 
2-loop reduced coupling. The requirement Bf < translates into the condition 



3R 3 - 12R 2 + 
which yields (for N 2 > 1/13), 



13 + w 



R - 3 < 0, (50) 



R& 2.314. (51) 

Note that this line is below the top of the window, and thus in the upper part of 
the conformal window the dual coupling has a space-like singularity. 



3.3 Summary 

To conclude this part, let us summarize the differences between QCD and SQCD 
with respect to the analyticity structure of the coupling in comparison with the 
boundaries of the conformal window (fig. 1). 

In QCD, the region of a causal 2-loop coupling covers the entire conformal win- 
dow (supposing the lower boundary is determined by super convergence: r = 3.25). 
As r is reduced further (below r ~ 3.17), the 2-loop coupling develops a couple of 
Landau branch points at complex Q 2 values. At even lower r, below r ~ 2.62, a 
Landau branch point appears on the space-like axis. 

Studying higher loop effects we showed that the 3-loop term is important in the 
lower part of the conformal window, and so the 3-loop coupling should be referred 
to as a zeroth order approximation in the infrared. The next observation is that the 
3-loop coefficient is negative in the conformal window both in MS and in all the 
physical effective charges for which the 3-loop coefficient has been computed. This 
means that the 3-loop coupling is causal at least where the 2-loop coupling is, i.e. in 
the entire conformal window, and in many cases, depending on /3 2 , also somewhat 
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below this region into the upper part of the confining phase. The 3-loop solution is 
reliable according to the usual perturbative justification: the 4- loop term in the f3 
function, at least in MS , is small enough not to affect the 3-loop solution. 

In SQCD, the region of a causal 2-loop coupling 1.8377 ^ R < 3 does not 
cover the lower part of the conformal window (the lower boundary is at R = 1.5). 
Below R ~ 1.8377 the 2-loop coupling develops a couple of Landau branch points at 
complex Q 2 , and below R ~ 1.5, i.e. below the conformal window (see eq. (38)) the 
2-loop coupling has a space-like Landau singularity. Studying higher orders we find 
that the 3-loop term is significant, and like in QCD it leads to a smaller coupling 
and to a larger causality region. But since the value of the coupling is still not small 
enough, and the 4-loop term is large, the 3-loop solution cannot be trusted. This 
means that the perturbative analysis in the electric theory is reliable only in the 
upper part of the conformal window. In the dual (magnetic) theory the reduced 
2-loop coupling is causal only in the region: 1.5 ^ it! ^ 1.8357. This coupling even 
has a space-like Landau singularity inside the window, for R > 2.314. 

Our main conclusions from this analysis are the following: 

(a) In QCD perturbation theory seems consistent in the infrared within the entire 
conformal window, and even somewhat below it. It then seems natural to as- 
sume that non-perturbative corrections are small, at least within the conformal 
window. 

(b) The previous assumption implies that in QCD the fields are, in some sense, 
weakly coupled even at the bottom of the window. This is contrary to SQCD 
where the electric fields are strongly coupled at the bottom of the window, one 
of the assumptions on which duality is based (see (d) below). We conclude 
that in QCD there is no dual description of the infrared in terms of some 
alternative degrees of freedom which are weakly coupled near the bottom of 
the window. 

(c) We found that the fixed-point in SQCD at the far-end of the conformal window 
cannot be explained in terms of the perturbative (3 function. 

(d) The regions where the electric and magnetic 2-loop couplings in SQCD arc 
causal do not overlap. Perturbation theory is never meaningful in the infrared 
in both the electric and magnetic descriptions of the same model. This is 
in accordance with the assumption on which duality is based that when the 
electric theory is weekly coupled, the magnetic is necessarily strongly coupled 
and vice-versa. 

(e) In SQCD perturbation theory signals its own inapplicability indicating that 
the coupling becomes strong within the window. This fits the same general 
philosophy on which the assumption in (a) is based: the strong coupling nature 
of the theory at the bottom of the conformal window should manifest itself 
already in perturbation theory. 
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4 Banks-Zaks expansion in SQCD vs. QCD 



In the previous sections we saw that in QCD perturbation theory yields a consistent 
description of the infrared physics even in the lower part of the conformal window: 
the coupling is causal and stable with respect to higher-loop corrections. On the 
other hand, in SQCD causality cannot be achieved at the perturbative level in the 
lower part of the conformal window. 

In order to examine the effect of higher order corrections we used an explicit so- 
lution of the equation (3{x) — in the MS scheme and in physical schemes in QCD, 
and in the DRED scheme in SQCD. Another natural way to study the value of the 
physical quantities in the infrared is the Banks-Zaks expansion, i.e. a power series 
solution to the equation j3(x) = 0, in terms of the distance from the top of the confor- 
mal window. In QCD, the expansion parameter is e = (11/2) — (Nf/N c ) = 6P /N C 
and the expansion has the form: 

x FP = z x e + z 2 e 2 + z 3 e 3 H (52) 

where Zi are independent of Nf. Since the coefficients of the (3 function are polyno- 
mials in Nf, it is possible to write them as follows. The 2-loop coefficient: 

c = 7r = -— + c lo ( 53 ) 

where a is proportional to e (and to /3 ) and Ci )0 is independent of Nf. The 3-loop 
coefficient: 

c 2 = = c 2 -l h c 2 , + c 2 ,ia + c 2)2 ao, (54) 

Po a o 

where c 2i j are independent on Nf, and so on. Then the leading terms in the Banks- 
Zaks expansion for a generic effective charge are [14, 15], 

x FP = a + (ci )0 + c 2 _i) al + ■■■ (55) 

We identify ao = z\t and note that z\ is the same for any effective-charge (or 
coupling) due to the universality of c. However, already z 2 depends on the effective- 
charge (or coupling) under consideration - according to eq. (55) it depends on the 
3-loop coefficient of the effective-charge (3 function. 

We stress that the ultimate justification of the presence of a fixed-point near 
the top of the conformal window, and thus of the very existence of the conformal 
window, is through this expansion [1, 14]. On the other hand, it is a priori not at 
all clear how far into the conformal window one can trust the expansion. We will 
be interested in particular in calculating the coupling and the critical exponent in 
QCD at the bottom of the conformal window and in estimating the reliability of this 
calculation. We will show that a calculation of this sort cannot be done in SQCD 
in the lower part of the conformal window. 
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4.1 Banks-Zaks expansion for the coupling 

4.1.1 Banks-Zaks expansion for the coupling in QCD 

As in the previous sections we start by considering the MS scheme. The advantage 
is that the coefficients of the (3 function are known up to 4-loop order [25]. This 
will enable us to compare the infrared limit obtained from the explicit solution of 
/3(x) = (fig. 3) which seems quite reliable at the 3-loop and 4-loop orders, to that 
of the Banks-Zaks partial-sums. A disadvantage of this scheme is that the coupling 
constant is not directly related to any measurable quantity. The dependence of 
the Banks-Zaks expansion on the effective charge or coupling under consideration, 
which was investigated in [11], first appears at the next-to- leading order term in the 
expansion - see eq. (55). This dependence becomes significant at the next-to- next- 
to- leading order level. 

According to [11] the next-to- next-to-leading order coefficient in the Banks-Zaks 
expansion for the MS coupling is rather large, making the corresponding term in 
the expansion comparable to the leading order terms already within the conformal 
window. Here we shall further analyze the expansion for the MS coupling explaining 
the source of the large next-to-next-to- leading coefficient. For physical effective 
charges this coefficient is smaller than in MS , hence the expansion is more reliable. 

The coefficients of the j3 function in the MS scheme are known up to 4-loop 
order [24, 25]. The three first Banks-Zaks coefficients in the expansion of xfjg (52) 
are then determined: 

16 N c 

Zl 3 25iV 2 - 11 1 ' 



16 N t 

%2 = = 



548 N c 4 - 1066 N c 2 + 231 



27 (25iV c 2 -ll)' 
32 N r J 



Z3 



243 (25iV c 2 -ll 



5 



with 



J = 52272 + (389235 - 1341648Cs)iV c + (-719758 + 3362832C 3 )iV c 4 
+ (-1148400C 3 - 1105385)iV c 6 + (990000C 3 + 730529)iV c 8 

Let us examine whether the Banks-Zaks expansion (52) is still reliable at the 
bottom of the conformal window. Table 1 summarizes the results for N c x^ (this 
normalization is used in order to consider both finite N c cases and the large N c 
limit) according to (52) and (56) at the lower boundary of the conformal window, 
namely at e = 11/2 - 13/4 = 9/4 for N c = 2, N c = 3 and N c — ► oo. The results 
are presented as a function of order in e: order e stands for the leading term in (52), 
order e 2 stands for the sum of the first two terms in (52), an so on. 
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order 


N c = 2 


iV c = 3 


N c — > oo 


e 


0.539 


0.505 


0.480 


e 2 


0.620 


0.601 


0.585 


e 3 


1.03 


0.933 


0.880 



Table 1: N c x^ in QCD at the bottom of the conformal window as a function of 
order in the Banks-Zaks expansion. 

Our first conclusion from table 1 is that there is no significant dependence on 
N c : there is no much difference between N c x^ for N c = 2 and for N c = oo. 

As mentioned above, the C(e 3 ) term at the bottom of the window is larger than 
the 0(e 2 ) term there. Note that it is also comparable to the leading 0(e) term. 
This clearly raises doubts concerning the reliability of the expansion. On the other 
hand, solving explicitly j3(x) = we found in Sec. 2 that the 4-loop fixed-point 
value is almost identical to the 3-loop one down to the bottom of the conformal 
window (fig. 3). This calls for a more detailed examination of the relation between 
the Banks-Zaks expansion and the explicit solution, which we conduct in the next 
section. 

4.1.2 The reliability of the Banks-Zaks expansion in QCD at the bottom 
of the conformal window 

The purpose of this section is to understand the reason for the large 0(e 3 ) term 
in the Banks-Zaks expansion in MS , and finally to estimate the reliability of the 
fixed-point value. The analysis we present is for the case N c — > oo, but the results 
for low N c are qualitatively the same. 

Let us compare first the numerical values obtained at the bottom of the window 
from the explicit solution vs. the corresponding partial sum in the Banks-Zaks 
expansion: 



order 




P(x) = 




e 


0.480 


2-loop 


2.1818 


e 2 


0.585 


3-loop 


0.7495 


e 3 


0.880 


4-loop 


0.7667 



Table 2: N c x^ in large N c QCD at the bottom of the conformal window as a 
function of order in the Banks-Zaks expansion and from an explicit solution of the 
equations j3{x) = for the truncated (3 function at each order. 

This comparison is shown also in fig. 8. We see that the two calculation proce- 
dures agree. Referring to the explicit solution as the best estimate at hand, we can 
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estimate the uncertainly in the value of the infrared coupling from the difference 
between the two calculation procedures. For the xfjg the uncertainty is no more 
than ±25%. 

Let now investigate the relation between the explicit solutions and the Banks- 
Zaks expansion. At 2-loop order, the functional form of the fixed-point value in the 
large N c limit is 

^ (0) = JVc (zi) = _J5l_. (57) 

At higher loop orders, the result is a more complicated function of e. At any order 
the explicit solution has a finite convergence radius in powers of e, and thus we 
expand it, and compare the expansion to the function itself. Such a comparison is 
shown in fig. 9 at the bottom of the conformal window, i.e. for e = 2.25. 

In the upper plot, corresponding to the 2-loop case, we see that the expansion 
in e converges very slowly to the explicit solution. This can be understood knowing 
that N c x(0) is a geometrical series in e (57) and that e at the bottom of the window 
is already quite close to the convergence radius which is e = 75/26 ~ 2.88, the point 
where c vanishes. Since we know from the comparison with the explicit solutions 
at higher orders that close to the bottom of the conformal window the 2-loop value 
for A^a^O) is unrealistically large§ we should not regard the slow convergence of the 
series in e corresponding to (57) as indicative of a problem of the Banks-Zaks series 
as a whole. It just means that higher orders are important. 

In the 3-loop case in fig. 9 (middle plot) the Banks-Zaks partial sum at order 
0(e 2 ) is much closer to the explicit solution and the convergence at higher orders in 
e is much accelerated as compared to the 2-loop case. 

In the 4-loop case in fig. 9 (lower plot) the partial sums of the e expansion diverge 
badly beyond the 0(e 3 ) term or so. The reason is that the convergence radius of 
the e series of the explicit solution is about e ~ 1, i.e. significantly smaller than 
e = 2.25 which corresponds to the bottom of the window and to fig. 9. This also 
explains why the 0(e 3 ) term in the Banks-Zaks series, which is fully determined at 
the 4-loop level, is larger than the 0(e 2 ) term. The explicit solution is a well defined 
function of e in the entire conformal window in all the cases considered. It turns 
out however that in the 4-loop case this function does not have a converging power 
expansion beyond e ~ 1. This fact is shown also in fig. 3: around e ~ 1 the series 
departs from the explicit solution itself. 

We note that for the available examples the e series that correspond to increasing 
loop-order solutions have an ever decreasing convergence radii: it is e ~ 2.88 in the 
2-loop case, e ~ 2.787 in the 3-loop case and e ~ 1 in the 4-loop case. This may be 
related to large order behavior of series: since the Banks-Zaks expansion is based on 
the factorially growing perturbative coefficients, it is natural to expect that it is also 
an asymptotic series with zero radius of convergence. Such a behavior will be avoided 
only if some systematic cancellation of the factorially growing ingredients occurs. 
If indeed the asymptotic nature of the Banks-Zaks series is reached at the order 

§ This is related to the discussion in Sec. 2 concerning the necessity to start from the 3-loop 
term in order to establish perturbative causality in the lower part of the window. 
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C(e 3 ) the best estimate of the fixed-point value from the expansion is obtained by 
truncating the series after the minimal term, in this case, the next-to-leading term. 

A comparison between the fixed-point value from the Banks-Zaks expansion and 
the explicit solution of (3(x) = can be also conducted in physical renormalization 
schemes. In the absence of full 4-loop perturbative coefficients, one cannot obtain 
an explicit solution at the 4-loop level. On the other hand, the 0(e 3 ) is calculable 
[14, 15, 11] and thus the next-to-next-to-leading order partial sum can be compared 
with the explicit solution of the 3-loop effective charge /3 function. Such a comparison 
was performed in [11] for the effective charge which is defined from the vacuum- 
polarization D-function. As shown in fig. 7 there, the two calculation methods 
nicely agree down to the bottom of the conformal window (iVj ~ 10 in the figure) 
and even below. 

As noted above, in physical renormalization schemes the Banks-Zaks coefficients 
(and in particular the next-to-next-to-leading coefficients) are smaller than in MS 
[11, 15], and so the expansion seems more reliable. For example, for N c = 3 we have 
[11]: 

x ™ = a + 1.14 ajj + 23.27 ajj + ••• (58) 

= a + 1.22 a 2 + 0.23 a 3 , + ••• 
xl v = a - 0.86 a 2 + 10.99a 3 + ••• 

where D stand for the effective charge defined from the vacuum polarization D- 
function and V stands for the one defined from the heavy quark potential. In fact, 
the Oiflo) coefficient in x pp is the largest amongst all the 0(ajj) coefficients for the 
effective charges considered in [11]'. 

We conclude that calculation of infrared quantities can be performed either as 
an explicit solution of the equation (3(x) =0 or by the Banks-Zaks expansion. Al- 
though the expansion probably has a zero convergence radius in general, and bad 
convergence properties already for the available 4-loop example (MS ) , it seems to 
give a reasonable estimate at the next-to-leading and the next-to-next-to-leading 
orders within the entire conformal window. Infrared quantities appear to be pertur- 
batively calculable in general even at the bottom of the conformal window. Note, 
however, that the accuracy is observable dependent. Some quantities, like the vac- 
uum polarization D-function, can be determined with high accuracy, whereas for 
others the accuracy is not as good: as mentioned above, the MS coupling can be 
determined within ±25% accuracy. 

4.1.3 Banks-Zaks expansion for the coupling in SQCD 

Let us now turn to the supersymmetric case and consider the Banks-Zaks expansion 
for the value of the DRED coupling at the fixed-point. The expansion parameter is 

^The result presented above for the O(ofy coefficient in x PP is different from the one in [11]. 
The latter was calculated based on a wrong 2-loop coefficient, which has now been corrected thanks 
to [26]. 
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5 = 3- R = 3- (N f /N c ): 



< RBD = Z x 8 + Z 2 5 2 + Z 3 5 3 + 0(5 4 ). (59) 

The coefficients of the (5 function up to 4-loop are taken from [27]. The resulting 
Banks-Zaks coefficients read: 

2 N c 

z ' = SivFT (60) 

7 1 JV, 

Zj2 — - 



3 (iV c 2 - 1) 

1 iV c [(17 + I8C3R 4 + (-25 + 18C 3 )iV c 2 + 8] 
3 ' 54 (iV c 2 -l) 3 

Table 3 summarizes the results for N c x^ ED , according to (59) and (60), at the 
bottom of the conformal window, i.e. at 5 = 3 — 3/2 = 3/2: 



order 


N c = 2 


N c = 3 


N c — > 00 


5 


1.33 


1.13 


1 


5 2 


2.33 


1.97 


1.75 


5* 


8.00 


5.38 


4.16 



Table 3: N c x^ ED in SQCD at the bottom of the conformal window as a function of 
order in the Banks-Zaks expansion. 

There is a clear contrast between the supersymmetric case of table 3 and the 
non-supersymmetric case of table I. Table 3 shows that the Banks-Zaks series for 
-^c^dred & t ^ ne bottom of the conformal window cannot be trusted at all, since the 
next-to-leading term is comparable to the leading one and the third order term is 
much larger than both. In addition, the value of the coupling itself (as much as it 
can be determined) is larger than in QCD. 

It is interesting to compare between the explicit solutions to the equations (5{x) = 
at increasing loop order (fig. 6), and the Banks-Zaks expansion. In the following 
table we show the values of the infrared coupling at the bottom of the window as 
determined by the two methods: 
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order 




ft(x) = 




5 


1 


2-loop 


oo 


P 


1.75 


3-loop 


4 


5 :i 


4.16 


4- loop 


no solution 



Table 4: N c x^ w in SQCD for N c — > oo at the bottom of the conformal window 
as a function of order in the Banks-Zaks expansion, and from the explicit solution 
of ft(x) = 0. 



It is clear from this table and from fig. 6 that the perturbative analysis fails 
to determine the infrared value of the coupling in the lower part of the conformal 
window. It thus seems, also from this point of view, that perturbation theory is 
inapplicable to describe the infrared physics there. 



4.1.4 Banks-Zaks expansion for the coupling in the magnetic theory 
(dual SQCD) 

In a similar manner we consider the Banks-Zaks series in the dual theory, where the 
expansion parameter is 6 a — R — (3/2) H, 

*duai = Z% + Z% + Z d 3 5j + 0(Si). (61) 

The coefficients can be calculated either from the reduced ft function (46), or directly 
from the coupled ft function (40), assuming both infrared couplings are vanishingly 
small. Since the ft function in the magnetic theory is known at present only up to 
the next-to-leading order term, only the leading order coefficient in the Banks-Zaks 
expansion can be calculated. The result is: 

112 

1 3 iV c 2 -4 y ! 

The infrared value of the Yukawa coupling is given by 

a fp = YW e Sd + ° { ®- (63) 

Let us now examine the magnitude of the infrared coupling in the magnetic 
theory at the top of the conformal window (having only one term, we cannot in- 
vestigate the behavior of the series as we did for the electric theory and for the 
non-supersymmetric case). Using the leading term in (61) with (62) and 5d = 3/2 
we find that for N c = 3, x FP ~ 168/5 = 33.6, and for N c = oo, A^ c x FP ~ 56. In both 
cases, it is clear that the coupling is much too large to be perturbative (which also 

"Note that both expansion parameters S and Sd are chosen to be positive inside the conformal 
window. 
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implies that these values are meaningless). The conclusion is that the fixed-point 
of the dual theory cannot be described by perturbation theory at the far-end of the 
window. 

An interesting unrelated observation is that for N c = 2, the Banks-Zaks expan- 
sion is completely ill-defined due to the pole at N% = 4 in (62). In the absence of 
the Banks-Zaks expansion it seems hard to establish the existence of a fixed-point. 
In fact, as we explain below, the problem is specific to the point around which the 
expansion is done, and therefore it may not imply anything special for the rest of 
the conformal window for N c = 2**. The original theory in this case (at the bottom 
of the conformal window) is an SU(2) gauge theory with Nf = 3. The implied dual 
theory has a color group of Nf — N c — 1, which means that there are no gluons. 
Mathematically, this appears as an ill-defined expansion since the point where the 
next-to-leading coefficient of the /3 function Bf vanishes (see eq. (47)) coincides 
with the point where B$ vanishes^, and thus the ratio Bfi/Bf which is usually used 
to define the expansion parameter Sd is not arbitrarily small near the point B$ = 
but is finite there. 



4.2 Banks-Zaks expansion for the critical exponent 

The critical exponent 7 has a special status since it is a universal quantity [32] : it de- 
termines the rate at which any perturbative coupling or effective-charge approaches 
its infrared limit^. In addition, discussing the analyticity structure of the coupling 
we found that the value of 7 is indicative of a causal coupling. Thus it is interesting 
to study the Banks-Zaks expansion and its break-down for this particular quantity. 

Let us start with a brief review of the definition and the basic properties of 7*. 
The critical exponent is defined as the derivative of the (5 function, 

P(x) = -P x 2 (l + cx + c 2 x 2 + ■■■) (64) 

at the fixed-point: 



d(3(x) 

7 = 



dx 



-/3 x FP 2 + 3cx FP + 4c 2 (ajpp) 2 H (65) 



X=Ipp 



from which eq. (29) follows. 

As already mentioned 7 is universal, i.e. independent of the renormalization 
scheme. To be precise, this statement is true so long as the transformations relating 
the different schemes are non-singular (see ref. [32, 33] and appendix B in ref. [15] 
and references therein). 

**The authors are in debt to D. Ansclmi for explaining this point. 

ttpor any N c > 2, Bf becomes negative already at lower Nf/N c , before Bq vanishes. 

in QCD was discussed in various papers; see for instance [14, 15, 33, 11]. 
*The notation is again that of QCD, but the same equations are relevant in SQCD, with the 
replacement of e by 5, fli by B i: c by d, and so on. 
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The Banks-Zaks expansion for 7 can be calculated using (65) together with the 
Banks-Zaks series for x FP , yielding a Banks-Zaks series of the form, 

7 = 9l e 2 + g 2 e 3 + + ■■■ (66) 

Note that contrary to a generic effective charge, the expansion for 7 begins with an 
e 2 term. A further difference is that the coefficients of (66) have an additional factor 
of N c , as compared to those of (52). 

It was shown in [14] that the coefficients ^ are universal, i.e. they are the same 
for any effective-charge x. This is in agreement with what is expected on general 
grounds, since 7 itself is independent of the renormalization scheme in which the (3 
function is defined, and the expansion parameter e is a well defined physical quantity. 

An additional interesting observation [14] is that the first two terms in the Banks- 
Zaks expansion for 7 are determined from the 2-loop f3 function: 

7 = gi e 2 + g 2 e 3 + • • • = {g x j z\) [a 2 + c h0 a 3 + • • •] (67) 

where a = z±e and ci j0 are defined in (53). Since g 2 is fixed by the 2-loop (5 function 
which is the leading order in which the Banks-Zaks fixed-point can be discussed, it 
makes sense to regard the first two orders g\e 2 + g 2 e 3 together as the leading term. 
We shall see below that in both QCD and SQCD g 2 e 3 is comparable to g±e 2 for values 
of e such that the expansion for the coupling is still reliable*. However, according 
to the explanation above this should not be regarded as an indication of the break 
down of the series - it is the magnitude of the next term g 3 e 4 , that depends also 
on the 3-loop and 4-loop coefficients of the (3 function, which must be examined in 
order to assess the reliability the expansion. 

4.2.1 The critical exponent in QCD 

Again, we start with QCD where the coefficients of the Banks-Zaks series for 7 in 
(66) are*: 

8 iV 2 

91 = 9 25ivf3II (68) 
16 iV c 2 (13iV 2 -3) 



92 = 

9s = 



27 (25iV 2 -ll) 2 

8N 2 H 
243 (25 iV 2 - ll) 5 



with 



H = 3993 + (571516 - 894432 ( 3 ) N 2 + (-1599316 + 2241888 ( 3 ) N*+ 
(-765600 Cs + 865400) iV c 6 + (660000 Cs - 366782) iV c 8 



tin QCD it is the case in all the physical renormalization schemes that where examined in [11], 
since always C2,-i < 0. Thus it turns out that the next-to-leading coefficient in (55) is smaller in 
absolute value than the one in (67). 

*The third order coefficient of the Banks-Zaks expansion for 7 has been calculated for the first 
time in [15]. 
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Our aim is to see whether 7 can be calculated from this expansion even at the 
bottom of the conformal window, and then with what accuracy. Fig. 5 (upper plot) 
shows, in addition to the results of the explicit calculation in various schemes, the 
following Banks-Zaks partial sums: gie 2 , gie 2 + g 2 e 3 , and g±e 2 + g 2 e 3 + ^e 4 as a 
function of e. The next-to-leading term is relatively large, but as explained in the 
previous section this should not be taken as an indication for the breakdown of the 
expansion. The relevant observation is that the next-to-next-to- leading term is just 
a small correction. At this level the Banks-Zaks series for 7 seems reliable. 

The comparison between the explicit calculation based on a truncated /3 function 
and the Banks-Zaks partial sums, shown in fig. 5 (upper plot) raises again the 
question of the relation between the two calculation procedures, especially in the 
MS scheme. 

The table below summarizes the numerical values obtained at the bottom of the 
conformal window (like the plot, the numbers correspond to N c — > 00 but the 
results at low N c are similar). 



order 




P(x) = 


MS 


D 


Bj 


Fi 


V 


e 2 


0.180 














e 3 


0.320 


2-loop 


0.818 


0.818 


0.818 


0.818 


0.818 






3-loop 


0.466 


0.330 


0.340 


0.337 


0.413 


e 4 


0.284 


4-loop 


0.463 











Table 5: 7 in large N c QCD at the bottom of the conformal window as a function of 
order in the Banks-Zaks expansion and as an explicit calculation from the truncated 
(5 function in MS and in various physical schemes: D - vacuum polarization D- 
function, Bj and F\ - polarized and non-polarized Bjorken sum-rules and V - heavy 
quark potential. The gaps in the table are due to the fact that the next-to-leading 
term 0(e 3 ) of the Banks-Zaks series depends only on the 2-loop f3 function, while 
the next-to-next-to-leading term 0(e 4 ) is determined by the 4-loop f3 function [14, 
15, 11]. 

Considered separarately, both the Banks-Zaks expansion and the explicit cal- 
culation in MS seem reliable. Still the disagreement between them is about 40%. 
In order to understand better the source of this discrepancy we compare in fig. 10 
the explicit results for 7 with the partial sums in the e expansion of these results, 
at the bottom of the window. In the 2-loop (upper plot) and 3-loop (middle plot) 
cases the e series converges to the value for 7, while in the 4-loop case, the series 
diverges since its convergence radius is smaller than the value of e at the bottom of 
the window, e = 2.25. 

The comparison in fig. 10 suggests that the explicit calculation (right column in 
the table) is equivalent to some resummation of higher order terms in e, and explains 
the disagreement between the two calculation procedures. Such a resummation is 
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necessarily scheme dependent since it should reflect the spread between the differ- 
ent schemes when using a truncated (3 function. Finally, at the available order in 
perturbation theory we can determine the critical exponent to be 7 = 0.4 ±0.1. 



4.2.2 The critical exponent in SQCD 

In the SQCD case, in the original (electric) theory, the expansion for the critical 
exponent is 

7 = G^ 2 + G 2 5 3 + G 3 5 4 + ■■■ (69) 

and the coefficients Gj are: 

1 N 2 

1 N 2 (2N 2 - 1) 
18 (AT? - l) 2 

1 N 2 [(-1 + 18C 3 )iV c 4 + (2 + 18C 3 )iV c 2 - 5] 
216 (.V-' 

The numerical values of 7 as calculated from the partial sums Gi5 2 , G\5 2 + G 2 8 3 , 
and G\5 2 + G 2 5 3 + G^5 4 is shown in fig. 5 (lower plot) as a function of 5 within the 
conformal window. It is clear from the plot that the expansion is useless at the 
bottom of the window since G3S 4 is comparable to G\5 2 and to G 2 5 3 . 



G 1 = 
G 2 = 
G 3 = 



4.2.3 The critical exponent in the magnetic theory (dual SQCD) 

Finally we consider the Banks-Zaks expansion for the critical exponent in the dual 
SQCD theory, 

1 = G d 1 5 2 + G d 2 5 3 + G d 3 5 4 d + --- (71) 

There are two ways to calculate this quantity, one, which has been used in [35], is 
based directly on the coupled renormalization group equations (40) and the other 
is based on the reduced equation (46). We show that both methods give the same 
Banks-Zaks expansion. 

Calculating 7 in the magnetic theory directly from the coupled renormalizat ion- 
group equations (40) is more involved, since there are two couplings. As mentioned 
above, a similar calculation was performed in [35]. The latter ref. presents a calcu- 
lation of the leading-order term in the expansion, but in fact, as we shall see, the 
2-loop gauge (3 function together with the one-loop Yukawa (3 function fixes also the 
next-to-leading order term, just like in QCD and in the SQCD electric theory. 

Let us briefly describe the method and then give the results. The generalization 
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of 7 to a two coupling theory is the following matrix: 



/ d(5 d x d(5j \ 
dx d\ 



dpi 



dpf 



V dx dX I 



(72) 



The next step is to diagonalize the matrix. This yields two eigenvalues: 71 and 72. 
Therefore a physical quantity behaves in the infrared like 



+ K 2 



(QVAL) 



72 



(73) 



and then asymptotically only the minimal eigenvalue is important. Thus we conclude 
that 7 = min {7!, 72}. 

Taking the derivatives of the coupled P functions (40) at the fixed-point we find 
the matrix elements of (72): 



dPl 



dx 
dPi 



392 JV 



T2 



1120 iV 2 (4 + 13 iV 2 



9 (TV 2 - 4)= 



dX 

dpi 



= -784 



Nt 



dx 

dpi 



(N? - 4) 

m 2 -4) 

3 N? 



-A 



448 iV c 4 (-36 + 65iV c 2 ) 3 
— j 



(iV 2 - 4)= 



+ 



8 (13iV 2 + 28) 
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iV 2 



<8 + 



dX 

The eigenvalues are 

71 = 

72 = 



28 8 (79 TV 2 + 76) tf 



28 . 

y 5d + 9 



820 ^ 2 + 76 , 2 + ^3) 



iV 2 -4 



(74) 



56 



3 N} 



N 2 

4 d 



64 (17iV 2 + 2)N 2 C 
9 (iV 2 - 4) 2 



The two eigenvalues are positive reflecting the infrared stability of the fixed point. 
The smaller eigenvalue is 7 = 72 . 

We note that these eigenvalues do not agree with the leading order calculation 
in [35]. The reason^ is that ref. [35] uses the /3 function as it appears in [30] - see 
the comment concerning [30] below eq. (41). 

The second method to calculate 7 in the magnetic theory is to use the reduced 
renormalization group equation. Here we have a single coupling and thus the deriva- 
tive of the P function (46) at the fixed-point immediately yields the relevant 7. Per- 
forming this calculation we indeed find the same value as appears in (74) for 72. The 

§The authors arc in debt to D. Ansclmi for his help on this matter. 
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non-relevant perturbation corresponding to 71 does not even appear when using this 
procedure. 

As stressed above the critical exponent can be (and was) calculated up to the 
next to leading order term from the available coefficients in the magnetic theory 
(3 function. This conclusion is transparent in the reduction method: the 2-loop 
reduced (3 function is fully determined from the 2-loop gauge (3 function and the 

1- loop Yukawa (3 function. In particular it does not depend on the 2-loop Yukawa 
(3 function. Like in QCD and in the electric SQCD cases, it follows that 7 can be 
computed up to the next to leading order term. On the other hand, this conclusion 
is non-trivial when calculating 7 from the coupled (3 functions. It turns out that 
the next to leading order terms in the matrix elements of V, likewise in the larger 
eigenvalue 71, do depend on the (unknown) 2-loop terms of the Yukawa coupling (3 
function. On the other hand, the relevant eigenvalue 7 2 does not depend on these 

2- loop terms. 

We emphasize that the first method described above, i.e. to use the coupled (3 
function in order to define 7 as a matrix, and then take the minimal eigenvalue, 
is a completely general procedure. It is guaranteed that the second method that 
uses the reduced (3 function will also give the correct value of 7 once we choose 
the reduced solution that belongs to the infrared stable fixed-point, which we did. 
Note that there exists in this case another possible reduction that corresponds to 
A = [31, 34] and an unstable infrared fixed-point. Since A = means a zero 
superpotential, this reduction cannot correspond to the dual of an electric theory. 
Had we used this reduction instead of the relevant one we would have obtained 
different values for the infrared fixed-point and for 7. 

Contrary to QCD and the electric SQCD theory, higher order corrections to the 
coupled (3 functions of the magnetic SQCD theory are not known. As a result we 
cannot study the behavior and the break-down of the Banks-Zaks series beyond the 
next-to-leading order term in this case. 

4.3 Interpolating between the original theory and its dual 

The results of the previous sections indicate that in SQCD the electric theory is 
strongly coupled when the magnetic theory is weakly coupled and vise-versa. These 
results are in accordance with Seiberg's description of the conformal window. Since 
calculations are usually limited to the weak coupling regime, it is in general impossi- 
ble to compare between results obtained in the two theories. An interesting example 
of how the two theories can be compared for the infrared limit of a specific physical 
quantity, the total "hadronic" cross section ratio in e + e~ annihilation, R e + e -, is 
discussed in ref. [34]. 

We find it useful to express the relation between the infrared limit of physical 
quantities in terms of the relation between the corresponding effective charges. For 
instance, we can define an effective charge in the original (electric) theory, 

R e+e -{Q 2 ) = a + bx Re+e _(Q 2 ), (75) 
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where a and b are N c and Nf dependent constants, and another effective charge in 
the dual (magnetic) theory, 



R d +e _(Q 2 ) = a d + b d x d R + (Q 2 ). (76) 

e 

For a non-zero Q 2 the quantities in the two theories are not related, but in the 
infrared limit duality relates them. Close to the fixed-point we have: 

4. + ,_(Q 2 ) = 4, + ,_(o) -(£-)" 

According to duality, we have in the infrared limit, .Rg +e _(0) = i? e + e -(0), i.e. 

a + bx Re+e _ (0) = a d + b d x d Re+e _ (0) (78) 

and 

7 = Id- (79) 

The last equality was shown to follow from duality in [35] where 7 was identified as 
the anomalous dimension of the Konishi current, which is part of the superconformal 
algebra. A natural extension of duality would be to conjecture [36] that also the 
terms that describe the approach to the infrared fixed-point are the same in the two 
dual theories. In this case it makes sense to set the convention such that A eff = A effd . 

From Seiberg's description of the conformal window it seems reasonable to as- 
sume that taking into account the physical information from both the electric and 
magnetic descriptions together, we may be able to describe the infrared limit in 
the entire window. In the following we give an example how this can be achieved 
in practice by interpolating between the Banks-Zaks expansions in the electric and 
magnetic theories. 

The quantity we consider is the critical exponent 7 which was calculated in the 
original and the dual theories in the previous section. As explained above, the 
calculations in the two theories correspond to the same physical quantity. 

We choose to analyze the critical exponent in the large N c limit (a similar analysis 
is possible for any N c > 3). Using the expansions for 7 in electric theory (69) and 
in the magnetic theory 7 = 72 (74), we obtain in the large N c limit the following 
partial-sums, respectively: 

7 = i|3-Hf + i|3-^ + (-i <3 + _l_) [3-^ + ... (80) 
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Fig. 11 presents the functional form of 7, according to the above expansions. In 
the original theory we show Gi5 2 + G 2 5 3 and Gi5 2 + G 2 5 3 + G 3 5 4 , and in the dual 
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theory we show GfSj + G%5j. Each of the two expansions can be trusted just in 
some limited region around the expansion point. It is quite clear from this figure, 
especially if one compares the results obtained in both descriptions, that the per- 
turbative result cannot be extrapolated to the far-end of the conformal window. 
Thus, a straightforward comparison between the results obtained in the two dual 
description is impossible. On the other hand, it makes sense to interpolate between 
the two. 

We use the 2-point Pade approximants method to interpolate between the two 
series (80). The general idea is to construct a rational function which yields both 
the known series (80) in the original and dual theories when expanded in a Taylor 
series at R = 3 and R = 3/2, respectively. The calculation technique is explained 
in detail in chapter 8 in [21]. The resulting approximant is: 



This approximant is a rational polynomial of order [4/4]. Note that the numerator 
in (81) contains the double zero at both ends of the conformal window as implied 
by (80). In principle, there could be further possible [N/M] approximants, based on 
other rational polynomials. But, given the particular form of (80), other approxi- 
mants cannot be constructed at this order. 

The interpolating 7 PA of (81) is shown in fig. 11 together with the lines describing 
the partial-sums (80) which correspond to the Banks-Zaks expansions around R = 3 
and R = 3/2. We stress that in order to construct (81) we used nothing but 
the information contained in the coefficients of (80). A priori, a Pade pole could 
have appeared within the conformal window, which would probably mean that this 
interpolation technique is inappropriate. We find that such a pole does not appear. 
On the other hand, the convergence radius of the 8 and 5d expansions of (81) are 
rather small due to complex Pade poles. This should not be a surprise, as we expect 
the all order result to have zero radii of convergence in terms of S and Sd- 

Another interesting observation is that duality suggests an alternative way to 
define the expansion parameters for both the original and the dual theory. Referring 
to the original degrees of freedom, we defined R = Nf/N c , and then the expansion 
parameters were: in the original theory 5 = 3 — R, which is small at the top of the 
window and in the dual theory 5d = R — (3/2), which is small at the bottom of the 
window. If we instead start with the dual theory, then we consider the ratio 
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[23814 Cs + 10881 - (47628 Ca - 22734) R + (34398 Cs + 17685) R 2 

-(10584 Cs + 6216) R 3 + (1176 Cs + 880) R 4 } 



R = N } /N d c 



N f R 



(82) 



N f -N c R-l 



and define the following expansion parameters: 
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It is then possible to repeat the calculation of the Banks-Zaks coefficients in terms 
of S in the electric theory and Sd in the magnetic theory. An alternative way to 
calculate the expansion in terms of S and Sd would be simply to use the expansion 
of (80): substituting 5 and Sd in terms of S and Sd, according to (83) in (80) and 
Taylor expanding to the maximal order to which the coefficients are fixed (order S 4 
in the electric theory and order 5% in the magnetic theory) we find: 
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Note that the functional form and therefore the numerical values of the truncated 
expansions for 7 using the S and 5d variables is different (at any finite order) from 
those using the S and 5d variables. On the other hand, 7 is a physical quantity 
and therefore it cannot depend on such an arbitrary choice of expansion parameter. 
This is a disadvantage of using truncated series: they do not respect the invariance 
property of the full function j(Nf/N c ) to the choice of expansion parameter. 

Finally, we come back to the 2-point Pade Approximant (81) and ask how this 
function behaves under this change of expansion parameter. One can construct a 
[4/4] 2-point Pade approximant starting from the series (84) and compare it, as a 
function of Nf/N c to the [4/4] 2-point Pade approximant of (81). One would find 
that the two functions are identical. The reason is that the transformation relating 
R to R, and then also the transformations relating 5 and 5d to S and Sd, are all 
homographic transformations of the argument of the Pade approximant. It is then 
guaranteed by a mathematical theorem (see [21, 37]) for each expansion separately, 
that the diagonal Pade function is invariant. Note that this theorem holds in our 
example which is a diagonal [4/4] rational function, but it does not hold for non- 
diagonal rational polynomials. It is interesting to mention that the same invariance 
property of Pade approximants was shown to be significant in a different context 
[37], where Pade approximants are used to resum perturbative series in QCD. 

This invariance property of the diagonal 2-point Pade approximants with respect 
to the choice of expansion parameter, indicates that it is a good candidate to serve 
as an approximation to the function r y(Nf/N c ). 

In order to have a rough estimate of the accuracy of this approximation pro- 
cedure, we suggest the following: let us construct a higher approximant based on 
a guess for the higher order coefficients of the expansion. We choose to "add" a 
coefficient in the Banks-Zaks expansion in the magnetic theory, C3. Not knowing 
anything about higher order coefficients of the coupled (3 functions that are needed 
to determine this coefficient, we just use Gf = 0. We do not expect it to be a good 
guess (especially considering the presumably asymptotic nature of the expansion) 
but we use it just to check the sensitivity of the calculation. With a value of 
at hand, we can obtain a (non-diagonal!) [4/5] rational polynomial that reproduces 
the first three terms in each expansion: the first three correct coefficients in the 
electric theory, and the first two correct coefficients in the magnetic theory, with 
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the third coefficient put to zero. The resulting [4/5] approximant is also shown in 
fig. 11. We refer to the difference between the [4/4] approximant of (81) and the 
latter as a rough estimate of the error of the interpolation which is due to the lack 
of knowledge of higher order corrections. 

This exercise suggests that although the partial-sums can be considered to be 
a good numerical approximation to the physical quantity 7 only in some limited 
domain around the boundaries of the conformal window, an interpolating Pade 
approximant can be a good approximation for 7 in the entire conformal window. 

Unfortunately, in the absence of non-perturbative calculations, it is impossible 
to validate duality, nor to check our prediction for 7 in the interior of the conformal 
window. This calculation can surely be improved if higher-order terms in the (3 
function of either the electric or the magnetic theory will be available. It is of course 
possible that there exist non-perturbative corrections that limit in principle the ac- 
curacy of this calculation. One should note, however, that although the Banks-Zaks 
expansions we started with (80) are presumably asymptotic series, with zero radius 
of convergence, the resulting function (81) is well-defined in the entire conformal 
window. Thus combining the information from the two theories in this 2-point Pade 
method, is in fact also a way to resum the divergent series. 

If there were many known terms in the electric Banks-Zaks expansion one could 
use some resummation technique of this single series (such as Borel resummation or 
the application of an ordinary Pade approximant), hoping to get a vanishing value 
for 7 at the other side of the conformal window, thus adding a further consistency 
check to the duality conjecture. 

It is important to stress that if there was a phase transition somewhere within 
the conformal window, controlled by non-perturbative effects that are inaccessible 
from neither side of the conformal window, our interpolation procedure would not 
have been meaningful. However, the phase structure picture drawn in [16] suggests 
a completely smooth transition between the top of the conformal window where the 
electric theory is weakly coupled and the bottom of the window where the dual 
theory is weakly coupled. 

We emphasize that the interpolation for 7 serves here just as an example, and 
the method can be applied in general to any physical quantity that can be calculated 
in perturbation theory in both the original and the dual theories, such as the one in 
eq. (78). 

5 Conclusions 

The purpose of this paper is to understand the nature of the non-trivial infrared 
fixed-point that appears in asymptotically free theories such as QCD and SQCD 
if the number of light flavors is large enough. The main question we deal with is 
whether this fixed-point always originates within perturbation theory, or is it due 
to non-perturbative physics. 

Since more is known on the phase structure of SQCD, it is natural to discuss 
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this theory first. The argument for the presence of a fixed-point in SQCD in [16], 
just like in QCD [1], is a purely perturbative one based on the 2-loop /3 function. 
On the other hand, Seiberg's picture of the conformal window [16] assumes that the 
electric theory is strongly coupled near the lower boundary of the conformal window, 
and thus the fixed-point there is non-perturbative. Therefore, our first task dealing 
with the SQCD case was to confirm that the electric theory is strongly coupled at 
large distances near the bottom of the conformal window. Indeed we showed by 
considering the analyticity structure of the coupling and by examining the Banks- 
Zaks expansion, that the presence of the fixed-point at the far-end of the conformal 
window cannot be established in perturbation theory: it is a non-perturbative fixed- 
point. Our investigation further shows that the perturbative analysis signals its 
own inapplicability to describe the infrared already before the dual theory becomes 
weakly coupled. 

In order to describe the infrared physics in the entire conformal window it is 
useful to combine the information from perturbative calculations in both the electric 
and magnetic descriptions. We suggest to use an interpolation method - a 2-point 
Pade approximant - that combines the Banks-Zaks expansions around the two ends 
of the conformal window into a single formula. This method was demonstrated 
in Sec. 4.3 for the critical exponent, but in fact it is applicable to any physical 
quantity which is calculable in perturbation theory in both theories. The invariance 
property of the diagonal 2-point Pade approximant under the change of expansion 
parameter (a change which is motivated by the symmetrical realization of duality in 
the conformal window) makes it appropriate to describe , ~f(Nf/N c ). Still, this is of 
course not an exact calculation: it is quite clear that further corrections, both such 
that are accessible by perturbation theory in either the electric description or the 
magnetic one, and eventually such that are not, do exist. 

The lower boundary of the conformal window in the non-supersymmetric case is 
not as well established. If the superconvergence criterion is only a sufficient condition 
for confinement and not a necessary one, the conformal window may be narrower 
than shown in fig. 1 (upper plot). On the other hand, there are evidence from lattice 
simulations [10] indicating that the confining phase transition occurs at lower Nf. In 
spite of these contradicting findings, we refer to the lower boundary of the conformal 
window as the one implied by superconvergence. 

Considering the analyticity structure of the coupling and the Banks-Zaks expan- 
sion, we find that perturbation theory in QCD is self consistent even at the lower 
boundary of the conformal window. The most important evidence supporting this 
conclusion are: 

(a) the 2-loop and 3-loop couplings have a causal analyticity structure in the 
entire conformal window. At 3-loop order this is guaranteed provided f3 2 < 0, 
which holds in MS and in all physical schemes for which 3-loop coefficients 
are available. 

(b) thanks to the small 4-loop term in the MS scheme, the standard perturba- 
tive justification holds down to the infrared limit, and the 3-loop solution for 
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Xm§ (Q 2 ) can be trusted in the entire Q 2 plane. 

(c) calculation of infrared quantities from the Banks-Zaks expansion are in reason- 
able agreement with the explicit solutions even at the bottom of the conformal 
window - see fig. 8 here and fig. 7 in [11]. 

A crucial assumption we have made is that non-perturbative corrections are 
small within the conformal window so long as they are not implied by inconsistency 
of perturbation theory. This assumption makes sense provided no chiral symmetry 
breaking takes place within the conformal window (this is known to be the case in 
SQCD). Usually when chiral symmetry is broken, one expects quarks to develop 
dynamical masses. Then at momentum scales below these masses quarks will de- 
couple from the dynamics. This would clearly invalidate the perturbative analysis 
that takes into account Nf massless quarks. In particular decoupling could make 
effective charges non-monotonous with scale and thus invalidate the notion of a sin- 
gle valued /3 function. In the absence of chiral symmetry breaking, and in particular 
within the conformal window, it is natural to expect that effective charges evolve 
monotonously with scale. Thus the monotonicity property conjectured for specific 
observables by various suggested generalizations [39] of the c-theorem [38] could 
actually be a generic property of all effective charges within the conformal window. 

Consistency of perturbation theory together with the assumption that also non- 
perturbative effects are small implies that QCD is, in some sense, weakly coupled 
even at the bottom of the conformal window. The weak coupling nature of the 
QCD might be a disappointing message to anyone wishing to extend the notion 
of duality beyond the supersymmetric case. On the other hand, the fact that the 
infrared coupling is basically controlled by the perturbative (3 function opens up the 
possibility to analyze various physical quantities in the infrared near the confining 
phase transition. An example of such analysis is that of [6] which explains the phase 
transition at the bottom of the window through decoupling of quarks due to chiral 
symmetry breaking. 

The analyticity structure of the coupling seems to be indicative of the reliability 
of perturbation theory in the infrared in both QCD and SQCD. A priori one might 
suspect that not much could be gained by analyzing the singularity structure of the 
coupling which is scheme dependent. It turns out, however, that several character- 
istics of the j3 function are generic. One simple example, which is crucial for our 
analysis, is the negative sign of fa f° r physical schemes (and MS ) in QCD in the 
entire conformal window. It is not clear to us how general this property is. 

The most important universal infrared property of the (3 function is the critical 
exponent which turns out to be the key parameter in the condition for a causal 
analyticity structure. The condition < 7 < 1 is necessary for a causal coupling at 
any loop order and even beyond perturbation theory, while in practice at the 3-loop 
order it is usually also sufficient. 7 = corresponds to a free theory, as obtained at 
the top of the conformal window. In the supersymmetric case, duality implies that 7 
is the same in both theories, and thus 7 vanishes also at the bottom of the window, 
where the magnetic theory is free. In terms of the electric theory the vanishing of 
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7 can be most simply interpreted as corresponding to a smooth change of the (3 
function as Nf/N c crosses its critical value: the /3(x) curves continuously change 
from curves that cross the x axis at two near-by points to ones that do not cross the 
x axis at all. The critical curve corresponding to (5{x) at the bottom of the window 
just touches the x axis from below: 7 = implies a double zero. 

In QCD we trust perturbation theory down to the bottom of the window and find 
that 7 7^ there. On the other hand, it is clear that for some effective charges, such 
as the one associated with the potential between heavy quarks, there is no scale 
invariant behavior at large distances in the confining phase and hence the fixed- 
point must disappear as the lower boundary of the conformal window is crossed. 
The disappearance of the fixed-point with a non-vanishing 7 implies that the above 
picture of a smooth change of the /3 function cannot be realized. In this case, when 
the lower boundary of the window is crossed there is a discontinuous change in the 
form of the (5 function for x > x FP . This discontinuity, however, does not imply a 
jump in the coupling x(Q 2 ) at any finite Q 2 but only at the infrared limit itself. 
Finding that in QCD 7 can be perturbatively calculated and that it does not vanish 
at the bottom of the conformal window is a further indication of the absence of a 
dual description of the infrared in terms of some other weakly coupled fields. 
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Appendix — The analyticity structure of the NSVZ coupling 



In Sec. 3 we analyze the singularity structure of the 2-loop coupling in SQCD. 
Another natural choice for the SQCD f3 function, at the same level of approximation, 
is the NSVZ form with the matter field anomalous dimension calculated to first order 
[28, 29]. It is interesting to see whether the condition for causality of the coupling 
with this choice for the /3 function agrees with that of the 2-loop choice. 

The NSVZ (3 function at this order is: 



(3{x) = -B x 2 i-g (85) 



where B is given in (32), 



N f N 2 - 1 

E 1 = -1 -± (86) 

3N C - N f 2N C y ' 

and 

D = NJ2. (87) 

The (3 function in the NSVZ form yields, of course, the scheme invariant 2-loop /3 
function of (31) upon expansion, with C\ — D — E\ (see eq. (33)). 

Exact integration of the (3 function (85) yields a coupling constant that can 
be written explicitly as a function of the scale parameter t = \n(Q 2 /A 2 ) using the 
Lambert W function*: 

x(Q 2 ) = 4- 



G exp 



E x l + GW(z) 

(88) 

1 / - B 0j 



g r 1 + 1* 



where W(z) is defined by [22]: 

W(z) exp [W(z)\ = z. (89) 

and 

G=l-§^. (90) 

Next, one should specify the branch of the Lambert W function such that asymp- 
totic freedom will be obeyed (x should be real and positive at large (space-like) 
Q 2 > and approach zero as t — > oo). Since we are interested in the asymptot- 
ically free case where B > 0, we find that also E\ > and thus the sign of G 
determines both the sign of z and its magnitude in the ultraviolet: 



*This solution is very similar to the one presented in ref. [12] for the 2-loop and Pade improved 
3-loop coupling. 
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(a) if G > 0, z > and approaches infinity in the ultraviolet (z — > in the 
infrared). Then the relevant branch on the space-like axis is the principle 
branch Wq(z), which is a monotonically increasing positive function of z for 
z > 0; Wo(z) — ► oo as z — > oo thereby assuring asymptotic freedom. G > 
is obtained whenever E\ > D > 0. The singularity structure in the Q 2 plane 
can be analyzed following the lines of [12]. There are two possibilities (for 
G > 0): if B /(E 1 G) > 1, i.e. if B > E 1 — D, there is a pair of Landau 
branch points at complex Q 2 values. On the other hand if B /(EiG) < 1, i.e. 
Bo < Ei — D, there are no Landau singularities and the analyticity structure 
is consistent with causality. 

(b) if G < 0, z < and approaches zero in the ultraviolet (z — > — oo in the in- 
frared). In this case, the relevant branch is W-i(z) which turns to minus infin- 
ity for z — ► CT, thereby assuring asymptotic freedom. The singularity struc- 
ture is simple: there is a single Landau branch point at z = — 1/e corresponding 
to a certain Q 2 > on the space-like axis (t sing = {E 1 /Bq)[1 — G + G ln(— G)]). 

Note that contrary to the 2-loop or Pade improved 3-loop couplings in QCD, a 
simple pole from the denominator of (88) does not appear, since whenever W is 
real, its sign is the same as that of G, and thus it is guaranteed that 1 + GW > 0. 

We conclude that the condition for a causal coupling is G > and B /(E 1 G) < 1, 
where in fact the more restrictive condition is the second one. Thus, the NSVZ 
coupling is causal whenever 

B <E 1 - D. (91) 

The condition (91) coincides with the condition for causality of the 2-loop coupling 
Bq < —C\. Similarly, the condition to have a space-like Landau singularity in the 
NSVZ coupling, D > E±, coincides with the corresponding condition in the 2-loop 
coupling: C\ > 0. 
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Figure 1: The conformal window in QCD (upper plot) and SQCD (lower plot) is 
shown in the Nf — N c plane. In both plots the (3q = line, separating the infrared 
free phase from the ultraviolet asymptotically free phase is drawn as a continuous 
black line. This line is the upper boundary of the conformal window. The lower 
boundary of the conformal window as implied by superconvergence is drawn in gray. 
In the SQCD case, this last line is also the line below which the dual theory becomes 
infrared free. In both plots, the (black) dashed line shows the lower boundary of 
the region where the 2-loop coupling has a causal analyticity structure. Below this 
line and above the dot-dash line there are complex Landau singularities. Below 
the dot-dash line there is a space-like Landau branch point. In the lower plot, we 
also show in gray the dual lines which describe the analyticity structure of the dual 
coupling constant. 
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Figure 2: The compact domain in the complex coupling plane which corresponds 
to the entire complex Q 2 plane according to the 2-loop j3 function in N c = 3 QCD 
with Nf = 10 and Nf = 9.7, i.e. just above the minimal value of Nf required for 
causality of the 2-loop coupling (Nf ~ 9.683). 
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Figure 3: The infrared limit of the 2- loop coupling and the 3- loop and 4- loop MS 
couplings in large N c QCD as a function of e = (11/2) — (Nf/N c ). The dashed line 
corresponds to the 20th order Taylor expansion in e of the 4-loop coupling. The 
continuous vertical line represents the bottom of the conformal window implied by 
super convergence and the dashed vertical line shows the 2-loop causality boundary. 
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Figure 4: The relative magnitude of the four leading terms in the large N c QCD (3 
function in the MS renormalization scheme in the infrared limit, for various values 
of r = Nf/N c , from the top of the conformal window (r = 5.5) down to the bottom 
(r = 3.25). The normalization in each plot is such that the leading order term is 1. 
This means that the second is: (3ix/(3 , the third: /3 2 x 2 /A) and the fourth: (3 3 x 3 /(3 . 
The value of the coupling x is calculated as an explicit solution of the equation 
(3(x) = 0, with the 3- loop (3 function. Using the coupling as a zero of the 4-loop (3 
function does not change much the results. 
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Figure 5: The critical exponent 7 in large N c QCD (SQCD) is shown in the upper 
(lower) plot as a function of e = (11/2) - (N f /N c ) (5 = 3- (N f /N c )), according 
to an explicit calculation from the truncated 2-loop (3 function and the 3-loop and 
4-loop (3 functions (the loop order is indicated by the arrows) in the MS (DRED) 
renormalization scheme as well as according to the three available partial sums 
in the Banks-Zaks expansion: leading order - dotted line, next-to-leading order - 
dot-dashed line, and next-to-next-to-leading order - dashed line. In the QCD plot 
we show also the explicit calculation of 7 at the 3-loop order in various physical 
renormalization schemes defined from the vacuum polarization D-function (D), the 
polarized and non-polarized Bjorken sum-rules (Bj) and the heavy quark effective 
potential (V). The vertical arrow in the SQCD plot shows the point where the 
positive real zero of the 4-loop (3 function disappears. The vertical line is the bottom 
of the conformal window according to super convergence and the horizontal 7 = 1 
line is the necessary upper bound of 7 for a causal analyticity structure. 
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Figure 6: The infrared limit of the 2-loop coupling and the 3-loop and 4-loop DRED 
couplings in large N c SQCD as a function of 5 = 3—(Nf/N c ). The continuous vertical 
line represents the bottom of the conformal window implied by superconvergence 
(or duality) and the dashed vertical line shows the 2-loop causality boundary. 
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Figure 7: The relative magnitude of the four leading terms in the large N c SQCD 
P function in the DRED renormalization scheme in the infrared limit for various 
values of R = Nf/N C) from the top of the conformal window (R = 3) down to the 
bottom (R = 1.5). The normalization in each plot is such that the leading order 
term is 1. This means that the second is: B\x/Bq, the third: B2X 2 / Bq and the 
fourth: B 3 x 3 /B . The value of the coupling x is calculated as an explicit solution 
of the equation j3{x) = 0, with the 3-loop (3 function. 
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Figure 8: The infrared value of the coupling at the bottom of the conformal window 
(Nf/N c = 13/4) in large N c QCD, calculated as an explicit solution of the equation 
(3{x) =0 (circles), compared with the corresponding order partial sum in the Banks- 
Zaks expansion (crosses). The horizontal axis is the loop order of the calculation: 
2, 3 and 4-loop results in MS are shown. 
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Figure 9: The infrared value of the coupling at the bottom of the conformal window 
(Nf/N c = 13/4) in large N c QCD, calculated as an explicit solution of the equation 
(3(x) = 0. Three cases are shown in the three plots, from top to bottom: the 2- 
loop (5 function and the 3-loop and 4-loop j3 functions in the MS renormalization 
scheme. The horizontal line is the value of the infared coupling calculated from 
f3(x) = 0, and the symbols represent the partial sums in the expansion in powers 
of e = (11/2) — (Nf/N c ) of this solution. The cross symbols represent these partial 
sums that will not be altered by inclusion of higher order corrections to the (3 
function, i.e. they represent the Banks-Zaks partial sums. 
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Figure 10: The critical exponent 7 at the bottom of the conformal window (Nf/N c = 
13/4) in large N c QCD, calculated explicitly from the truncated (3 function. Three 
cases are shown in the three plots, from top to bottom: the 2-loop (3 function and the 
3-loop and 4-loop (3 functions in the MS renormalization scheme. The horizontal 
line is the value of 7 calculated from the truncated (3 function, and the symbols 
represent the partial sums in the expansion of 7 in powers of e = (11/2) — (Nf/N c ). 
The cross symbols represent these partial sums that will not be altered by inclusion 
of higher order corrections to the (3 function, i.e. they represent the Banks-Zaks 
partial sums which are renormalization scheme invariant. 
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Figure 11: The critical exponent 7 in large N c SQCD, is shown as a function of 
R = Nf/N c in the conformal window. In the original theory (expansion around 
R = 3) we show two Banks-Zaks partial sums, according to (80): second order partial 
sum - dot-dashed line; third order partial sum - dashed line. In the dual theory 
(expansion around R — 3/2) we show the second order partial sum as a dot-dashed 
line. The continuous line shows the [4/4] interpolating 2-point Pade Approximant 
of eq. (81). The dotted line shows the [4/5] approximant which is based on = 
and serves as a rough measure of the error for the [4/4] interpolating function. 
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